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Abstract 

This  effort  explored  the  impact  of  various  topographical  scales  present  within  irregular 
surface  roughness  on  a  turbulent  boundary  layer.  Low-order  representations  of  a  surface 
replicated  from  a  turbine  blade  damaged  by  deposition  of  foreign  materials  were  generated 
using  singular  value  decomposition  to  decompose  the  complex  topography  into  a  set  of 
topographical  basis  functions  of  decreasing  importance  to  the  original  (“full”)  surface.  The 
low-order  surface  models  were  then  formed  by  truncating  the  full  set  of  basis  functions  at 
the  first  5  and  16  modes  (containing  71%  and  95%  of  the  full  surface  content,  respectively), 
so  that  only  the  largest-scale  topographical  features  were  included  in  the  models  while  the 
finer-scale  surface  details  were  excluded.  Physical  replications  of  these  surfaces  were  rapid 
prototyped  and  particle-image  velocimetry  was  used  to  measure  the  flow  over  these  surfaces 
in  the  streainwise-  wall-normal  plane  and  in  a  streainwise-spanwise  plane  deep  within  the 
roughness  sublayer.  A  16-mode  model  of  the  full  surface  faithfully  reproduced  most  of  the 
characteristics  of  flow  over  the  full  surface  for  both  developing  and  developed  flow  conditions. 
However,  both  models  failed  to  reproduce  important  details  of  the  Reynolds-shear-s tress- 
producing  events  within  the  roughness  sublayer. 


1  Introduction 

Roughness  effects  can  play  a  crucial  role  in  a  variety  of  practical  engineering  systems,  from 
internal  flows  such  as  those  through  oil  and  gas  pipelines  to  external  flows  like  those  over  the 
surfaces  of  turbine  blades  and  heat  exchangers.  In  some  instances,  surface  roughness  occurs  in 
isolated  regions  of  a  flow  surface  meaning  that  the  flow  will  be  intermittently  perturbed  by  one 
or  more  step  changes  in  surface  condition  from  smooth  to  rough  and  vice-versa.  Such  transitions 
in  surface  quality  can  inhibit  the  flow  from  attaining  a  self-similar  state  (termed  developed  flow 
herein).  On  the  other  hand,  there  are  other  applications  for  which  surface  roughness  occurs  con¬ 
sistently  along  the  entire  length  of  a  flow  surface  of  interest,  meaning  that  its  development  may 
eventually  attain  self-similarity  since  the  flow  is  not  intermittently  perturbed  by  step  changes  in 
surface  quality.  Regardless  of  the  roughness  scenario,  it  is  of  significant  interest  to  understand 
the  impact  of  surface  roughness  for  improved  modeling,  prediction  and  eventually  control  of 
practical  flow  systems  in  the  presence  of  such  effects. 

Indeed,  many  studies  have  addressed  various  aspects  of  surface-roughness  effects  on  wall 
turbulence  under  both  developing  and  developed  conditions.  However,  despite  the  fact  that 
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the  surface  quality  in  engineering  systems  can  be  highly  irregular  due  to  multiple  damage 
mechanisms  and  can  vary  from  acrodynamically  smooth  prior  to  deployment  to  significantly 
roughened  over  time  due  to  cumulative  degradation,  most  of  these  efforts  employed  idealized 
roughness,  like  sand  grain,  woven  mesh  and  patterned  arrays  of  rough  elements.  While  such 
topographies  are  relatively  easy  to  realize  in  a  laboratory  setting,  their  ordered  character  and 
narrow  spectrum  of  scales  are  often  counter  to  the  surface  characteristics  of  practical  flow 
systems.  For  example,  turbine  blades  suffer  significant  surface  degradation  over  their  deployment 
lifetime  due  to  a  multitude  of  mechanisms,  including  deposition  of  foreign  materials,  pitting 
and  spallation.  Close  inspection  of  aged  turbine  blades  reveals  surface  degradation  marked  by  a 
broad  range  of  topographical  scales  as  well  as  a  relatively  high  intermittency  of  occurrence  [1,  2]. 
Other  examples  of  irregular  roughness  include  accumulation  of  algae  and  barnacles  on  the 
exterior  surfaces  of  submarines  and  ships  [3]  as  well  as  cumulative  erosion  on  the  blades  of  wind 
turbines  operating  near  the  sea.  In  all  of  these  examples,  single-scale  roughness  arranged  in  ail 
ordered  manner  will  likely  not  be  representative  of  the  rich  topographies  encountered  in  practice. 
Such  differences  have  been  known  for  decades,  most  notably  discrepancies  in  friction  factor 
in  the  transitionally  rough  regime  between  Nikuradse’s  results  for  monodisperse  sand-grain 
roughness  [4]  and  Colebrook’s  relationship  based  on  “industrial”  roughness  containing  a  broad 
range  of  topographical  scales  [5,  6].  More  recently,  Bons  [2]  found  that  classical  scalings  for  skin 
friction  over  rough  walls,  derived  for  flows  over  idealized  roughness,  do  not  hold  for  some  flows 
over  scaled  turbine-blade  roughness.  As  such,  an  understanding  of  idealized  roughness  effects 
may  not  properly  extrapolate  to  the  more  practical  case  of  highly-irregular  surface  roughness. 

It  is  well-documented  that  an  abrupt  change  in  surface  conditions  from  smooth  to  rough 
leads  to  the  formation  of  an  internal  roughness  layer  that  grows  in  thickness  with  distance 
downstream  [7].  Studies  of  idealized  roughness  indicate  an  overshoot  in  the  wall  shear  stress 
just  downstream  of  this  step  change  in  surface  condition  [8,  9].  In  addition,  a  significant  en¬ 
hancement  in  the  production  of  turbulence  is  noted  in  the  immediate  vicinity  of  the  roughness 
which  yields  higher  values  of  both  Reynolds  normal  and  shear  stresses  within  the  internal  layer 
compared  to  the  upstream  smooth- wall  flow  [8,  9,  10].  In  contrast,  the  flow  outside  this  layer  re¬ 
mains  relatively  undisturbed.  Such  roughness  effects  diffuse  away  from  the  wall  with  increasing 
downstream  distance  from  the  step  change  in  surface  condition  until  the  internal  layer  even¬ 
tually  engulfs  the  entire  wall  normal  extent  of  the  boundary  layer  and  the  flow  approaches  a 
self-similar  (developed)  state.  Previous  studies  indicate  that  the  precise  growth  in  the  thickness 
of  the  internal  layer  is  tied  to  the  details  of  the  rough  surface  encountered  [9].  Finally,  roughness 
can  also  significantly  reduce  the  spatial  scales  of  the  flow*  within  the  internal  layer  for  developing 
flow  [8].  With  regard  to  more  irregular  roughness,  Wu  and  Christensen  [11]  studied  the  impact 
of  a  short  streamwise  fetch  of  turbine-blade  roughness  due  to  spallation  damage  and  marked  by 
a  broad  range  of  topographical  scales  on  incoming  fully-developed,  smooth-wall  turbulent  chan¬ 
nel  flow.  This  effort  revealed  both  the  formation  of  the  expected  internal  layer  upon  transition 
from  smooth-  to  rough- wall  conditions  as  well  as  enhanced  local  Reynolds  stresses  within  the 
internal  layer  of  this  developing  flow  due  to  large-scale  topographical  features. 

For  developed  flow,  wherein  the  internal  layer  has  growm  to  engulf  the  entire  boundary-layer 
thickness  and  the  flow  has  attained  a  self-similar  state,  it  is  wcll-accepted  that  roughness  governs 
the  character  of  turbulence  within  the  roughness  sublayer  (3  —  5 k  away  from  the  wall,  where 
k  is  a  measure  of  the  roughness  height).  Outside  the  roughness  sublayer,  however,  roughness 
may  not  have  a  direct  impact  on  the  flowr  as  many  studies  of  idealized  roughness  report  that 
the  turbulence  behaves  similarly  to  that  of  smooth-wall  flow  when  properly  scaled  [12,  13,  14, 
15,  16,  17,  18,  19,  20,  21].  These  observations  are  in  accordance  with  Townsend  [22]  who  first 
hypothesized  that  at  high  Reynolds  numbers  (Re),  the  turbulent  motions  in  the  outer  layer  are 
independent  of  surface  conditions  and  viscosity  except  for  their  role  in  setting  the  wall  shear 
stress,  rw  (and  hence  the  friction  velocity,  uT  =  [r^/p]1/2,  where  p  is  the  fluid  density)  and  the 
boundary-layer  thickness,  6 .  With  respect  to  rough-wall  flows,  this  hypothesis  implies  that  if  the 
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characteristic  roughness  height,  /c,  is  sufficiently  small  compared  to  6  then  the  direct  impact  of 
roughness  is  eonfined  within  the  roughness  sublayer.  Under  such  conditions,  the  turbulenee  in 
the  outer  layer  is  only  indireetly  influenced  by  roughness  through  its  role  in  determining  uT  and 
6  [23].  Previous  studies  [24,  19]  indicate  that  this  seale  separation  as  represented  by  S/k  must 
exeeed  40  50  for  outer-layer  similarity  to  exist  (where  k  is  either  taken  to  be  the  characteristic 
roughness  height  or  the  equivalent  sand-grain  height,  ks). 

Outer-layer  similarity  has  also  been  reported  for  turbulent  flow  over  more  praetieal  rough¬ 
ness.  For  example,  Allen  et  al.  [25]  studied  turbulent  pipe  flow  in  the  presence  of  a  honed  surfaee 
akin  to  the  industrial- type  roughness  of  Colebrook  [6]  and  reported  smooth-  and  rough-wall 
mean  velocity  defect  profiles,  streamwise  turbulence  intensity  profiles,  and  streamwise  velocity 
speetra  that  eollapsed  in  the  outer  layer  in  accordance  with  Townsend’s  wall  similarity  hypoth¬ 
esis.  Similarly,  Wu  and  Christensen  [26]  reported  that  the  turbulenee  statistics  outside  the 
roughness  sublayer  remain  unaffected  by  roughness  replieated  from  a  turbine  blade  damaged 
by  deposition  materials  (compared  to  smooth-wall  flowr)  when  one  accounts  for  the  increased 
drag  at  the  surfaee  when  sealing  the  statistics  (using  uT  and  5).  This  similarity  was  also  found 
to  extend  to  the  average  spatial  structure  of  the  flow  through  comparison  of  two-point  velocity 
correlation  coefficients  outside  the  roughness  sublayer.  Thus,  Townsend’s  hypothesis  provides  a 
simple  means  of  predieting  outer-layer  behavior  for  rough-wall  flows  based  simply  on  knowledge 
of  uT  and  5 ,  though  previous  studies  indicate  that  it  may  not  be  a  universal  characteristic  of  all 
developed  rough-wall  flows  [27,  28,  29,  30].  Further,  while  the  outer  layer  may  not  be  direetly 
influenced  by  the  details  of  the  roughness  topography  in  question,  such  details  will  undoubtedly 
have  a  defining  impaet  on  the  loeal  flow  behavior  within  the  roughness  sublayer.  Thus,  given 
the  erueial  importance  topographical  details  ean  play  in  both  developing  and  developed  flows, 
and  the  fact  that  idealized  roughness  characterizations  generally  embody  a  rather  restricted  dis¬ 
tribution  of  topographical  seales,  the  effect  of  idealized  roughness  conditions  upon  wall-bounded 
turbulence  may  be  insufficient  for  successful  modeling  and/or  eontrol  of  practical  flows  in  the 
presence  of  irregular  roughness. 

Given  the  topographical  complexity  of  realistic  roughness,  the  relative  impact  that  each 
topographical  scale  of  an  irregular  surface  has  on  the  flow  is  eertainly  of  interest.  That  is,  are 
flows  over  irregular  roughness  predominantly  governed  by  the  impact  of  the  largest  roughness 
seales  or  do  the  finer  surfaee  features  contribute  in  a  meaningful  way?  Colebrook  and  White  [5] 
reeognized  the  importance  of  this  issue  in  their  studies  of  industrial  roughness  and  reported 
an  enhancement  in  pipe-flow  frietion  factor  with  the  addition  of  sandgrain  roughness  to  larger 
roughness  protrusions  compared  to  flow  over  the  larger  protrusions  alone.  Sehultz  and  Flack  [31] 
compared  the  flow  over  uniform  spheres  and  the  same  topography  with  the  addition  of  finer- 
seale  sandgrain  roughness.  They  reported  good  agreement  between  the  Reynolds  stress  profiles 
for  flow  over  both  surfaces  throughout  the  boundary  layer,  indicating  little  effect  of  finer-scale 
roughness  on  the  turbulenee.  However,  as  with  the  experiments  of  Colebrook  and  White  [5],  it 
is  not  elear  how  this  study  of  two  roughness  scales  of  substantially  different  size  translates  to 
the  ease  of  praetieal  roughness  which  is  marked  by  a  broad  speetrum  of  topographical  scales. 
With  respect  to  more  irregular  roughness,  Itoh  et  al.  [32]  measured  turbulent  flow  over  the  fur 
surface  of  a  seal  which  exhibited  a  riblet-like  character,  though  both  the  amplitude  and  wave¬ 
length  of  these  topographical  features  varied  significantly  in  space.  Itoh  et  al.  [32]  also  made 
measurements  of  turbulent  flow  over  a  model  of  the  seal  fur  that  consisted  of  ordered  riblets 
manufactured  with  the  dominant  amplitude  and  wavelength  of  the  seal  fur.  Comparison  of  these 
results  indicated  a  lack  of  consistency  in  the  turbulenee  statistics  for  flowr  over  the  real  seal  sur¬ 
faee  and  the  model.  These  differences  highlight  the  importance  of  topographical  eharaeteristies 
beyond  simply  the  dominant  amplitude  and  wavelength  of  the  real  surfaee  in  determining  its 
impaet  on  the  flow.  Finally,  Johnson  and  Christensen  [33]  considered  the  development  of  low- 
order  topographical  models  of  roughness  replicated  from  a  turbine  blade  that  contained  deep 
reeesses  of  varying  size  due  to  spallation  damage.  Model  topographies  w^ere  developed  using 
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singular  value  decomposition  (SVD)  and  short  fetches  of  the  models  were  fabricated  and  tested 
in  turbulent  channel  flow  in  which  the  upstream  smooth- wall  flow  was  fully-developed.  Under 
these  developing,  internal  flow  conditions,  it  was  found  that  a  model  containing  only  the  larger- 
and  intermediate  scales  of  the  topography  (10%  of  the  total  modes)  adequately  reproduced  the 
single-point  statistics  of  flow  over  the  full  surface  within  the  internal  layer  formed  by  the  abrupt 
transition  from  smooth  to  rough  conditions. 

The  present  contribution  builds  upon  the  initial  efforts  of  Johnson  and  Christensen  [33]  by 
considering  the  development  of  low-order  topographical  models  of  a  different  highly-irregular 
surface  roughness  replicated  from  a  turbine  blade  damaged  by  deposition  of  foreign  materials. 
The  ability  of  these  models  to  reproduce  the  characteristics  of  flow  over  the  original  roughness 
was  then  assessed  under  both  developing  and  developed  flow  conditions  in  a  turbulent  boundary 
layer.  Singular  value  decomposition  was  used  to  decompose  the  highly-inhomogeneous  surface 
topography  into  a  set  of  basis  functions  of  decreasing  contribution  to  the  overall  topography. 
Only  the  most  dominant  of  these  basis  functions  were  used  to  reconstruct  the  surface  topography, 
meaning  that  a  substantial  fraction  of  the  larger-scale  surface  features  were  included  in  the  low- 
order  models  while  the  finer  topographical  details  are  neglected.  Short  and  long  stream  wise 
fetches  of  these  low-order  representations  were  then  fabricated  and  tested  in  a  zero- pressure- 
gradient  turbulent  boundary  layer  to  assess  how  well  they  reproduce  the  flow  modifications 
generated  by  the  full  surface  topography  under  both  developing-  and  developed-flow  conditions. 
Such  comparisons  are  meant  to  reveal  the  relative  importance  of  finer-scale  roughness  features 
compared  to  the  most  dominant  roughness  scales  in  the  context  of  highly  irregular  roughness 
in  an  external  flow  arrangement. 

2  Experiments 

Particle-image  velocimetry  (PIV)  in  the  streamwise-wall-normal  (x  —  y )  plane  and  stereo  PIV 
in  a  streamwise- spanwise  plane  deep  within  the  roughness  sublayer  at  y  =  0.047J  were  used 
to  study  the  impact  of  low-order  models  of  surface  roughness  replicated  from  a  turbine  blade 
damaged  by  deposition  of  foreign  materials  on  a  zero-pressure-gradient  turbulent  boundary  layer 
at  Re#  =  Ue6/v  «  15000  for  developing  flow  and  13000  for  developed  flow,  where  Uc  is  the 
free-stream  velocity,  9  is  the  momentum  thickness  and  v  is  the  kinematic  viscosity  of  the  fluid. 
All  measurements  were  made  in  an  Eiffel-type,  open  circuit,  boundary-layer  wind  tunnel  with 
a  documented  turbulence  intensity  of  0.16%  in  the  free  stream  [34].  The  boundary- layer  plate 
utilized  in  all  of  the  experiments  consists  of  two  3  m-long  by  1  m-wide  flat  plates  suspended 
above  the  bottom  wall  of  the  tunnel  that  are  smoothly  joined  at  the  streamwise  center  of  the 
test  section. 

2.1  Roughness 

The  rough-wall  conditions  considered  herein  were  derived  from  one  of  the  surfaces  characterized 
by  Bons  et  ai  [1]  and  subsequently  used  by  Bons  [2]  to  study  bulk  skin  friction  and  heat  transfer 
characteristics  over  turbine-blade  roughness  (surface  4  in  Bons  [2]).  It  should  be  noted  that  the 
original  profilometry  measurements  of  this  damaged  turbine-blade  surface  by  Bons  et  ai  [1] 
yielded  roughness  heights  on  the  order  of  tens  to  hundreds  of  microns.  Therefore,  in  order  to 
generate  fully-rough  conditions  for  the  relatively  thick  boundary  layers  generated  by  the  flow 
facility  employed  (J  ~  100  mm)  at  the  Re#  considered  herein,  the  original  profilometry  informa¬ 
tion  was  scaled  up  in  all  three  dimensions  to  yield  a  topographical  condition  with  k  =  4.25  mm 
(Following  Bons  [2],  the  characteristic  roughness  height,  k,  is  taken  to  be  the  average  peak-to- 
valley  height.).  The  topographical  features  of  this  surface,  shown  in  figure  1(a),  are  elliptical 
in  shape,  are  generally  aligned  in  the  streamwise  direction  and  are  attributable  to  cumulative 
deposition  of  foreign  materials  on  the  blade  surface.  However,  a  broad  range  of  topographical 
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Figure  1:  (a)  Topographical  map  of  the  turbine-blade  roughness  under  consideration  in  the  present 
effort  (the  “full”  surface).  (b,c)  Low-order  representations  of  the  full  surface  generated  with  the  first 
16  and  5  SVD  inodes,  respectively,  (d)  Streamwise  profiles  of  roughness  amplitude  for  the  full  surface 
and  the  two  low-order  representations  coincident  with  the  lines  in  (a) -(c)  which  demarcate  the  spanwise 
posit  ion  of  the  streamwise  wall-normal  PIV  measurement  plane  in  the  present  experiments.  (e,f)  Residual 
topographical  maps  for  the  16-  and  5-mode  topographical  models,  respectively. 


scales  is  also  clearly  evident  in  this  surface.  Of  particular  interest  ,  the  general  characteristics  of 
this  surface  are  representative  of  surface  roughness  encountered  in  many  practical  flows  of  in¬ 
terest  but  are  quite  distinct  from  the  “idealized”  roughness  typically  studied  in  the  laboratory, 
including  sand  grain,  mesh  and  ordered  arrays  of  elements.  Given  that  the  boundary  layers 
under  consideration  have  thicknesses  of  approximately  100  mm,  the  full-surface  condition  gives 
S/k  ~  24  —  30  (<5  is  taken  as  the  wall- normal  position  where  the  mean  streamwise  velocity  equals 
99%  of  the  free-stream  velocity  for  all  cases).  Further,  the  root-mean-square  (RMS)  roughness 
height,  krms,  for  the  full  surface  is  1.0  mm  while  its  skewness  and  flatness  are  0.16  and  2.27, 
respectively.  For  reference,  a  long  fetch  of  this  surface  was  studied  by  Wu  and  Christensen  [35] 
and  outer-layer  similarity  was  observed  in  all  single-  and  multi-point  statistics. 

The  low-order  reconstructions  of  this  surface  were  devised  using  SVD  of  the  fluctuating 
topography  about  the  mean  elevation.  This  methodology  is  akin  to  proper  orthogonal  decom¬ 
position  (POD)  as  it  provides  an  optimal  basis  for  describing  an  inhomogeneous  signal  of  interest 
(Only  the  details  of  SVD  pertinent  to  the  present  application  are  presented.  The  reader  is  di¬ 
rected  to  Chatterjee  [36],  for  example,  for  a  more  detailed  discussion).  To  begin,  the  fluctuating 
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surface  elevation  of  the  full  surface,  rjfu\\(x,z),  shown  in  figure  1(a)  can  be  decomposed  as 

L 

V  (x,z)  =  '%2aj<t>j(x,z),  (1) 

j= l 

where  the  afs  are  coefficients  of  the  expansion,  the  <f)j' s  represent  orthogonal  basis  functions 
and  L  is  the  total  number  of  basis  functions.  Using  this  decomposition,  one  can  approximate 
the  fluctuating  surface  elevation  of  the  full  surface  as 

A/ 

Vfu\\(x,z)  ^  ?7jow(x,  2)  ^  ^  2),  (2) 

i- 1 

where  Tj is  referred  to  as  the  low-order  (or  reduced-order)  representation  of  7/ fun  based  on  only 
the  first  M  (<  L)  modes  of  the  decomposition.  Therefore,  some  detail  of  the  original  surface  is 
lost  in  this  reconstruction. 

The  basis  functions  and  coefficients  of  the  expansion  in  eqn.  (1)  are  determined  using  SVD 
One  can  write  the  topographical  information  contained  in  r/fun(x,2)  in  the  form  of  a  matrix  A 
and  this  matrix  can  be  decomposed  by  SVD  as 

A  =  UEVt,  (3) 

where  U  and  V  are  orthogonal  matrices  and  E  is  a  diagonal  matrix  containing  the  singular 
values  of  A,  an,  arranged  in  decreasing  value.  Relating  this  method  back  to  POD,  the  squares 
of  the  singular  values,  A j  =  ot2,  represent  the  eigenvalues  of  AAr  or  A1  A  while  the  columns 
of  U  and  V  are  related  to  the  eigenvectors  of  AA  r  or  A7  A,  respectively  [36].  Built-in  SVD 
functions  in  Matlab  are  used  to  determine  U,  £  and  V  which  yielded  a  total  of  L  =  383  modes 
for  the  full-surface  topography  of  figure  1(a). 

As  described  in  Johnson  and  Christensen  [33],  the  lowT-order  representations  of  the  surface 
elevation  can  be  realized  by  truncating  the  singular  values  in  £  at  mode  M  while  singular  values 
M+l  through  L  are  rendered  zero.  This  modified  singular- value  matrix,  £a/,  was  then  utilized 
to  compute  the  low-order  topographical  representation,  Am,  as 

Am  =  UEmVt,  (4) 

where  A  a /  embodies  77^  from  eqn.  (2).  Since  the  singular  values  are  arranged  in  descending 
order,  the  first  M  modes  represent  the  most  dominant  modes  of  the  surface  elevation  which, 
as  is  the  case  in  POD,  embody  the  larger  spatial  scales  of  the  original  topography.  Further,  as 
described  in  Johnson  and  Christensen  [33],  the  RMS  roughness  height  for  an  A/-order  model, 
k^p3,  can  be  related  to  feJJJf  as 

(klT)2  =  (5) 

l^j- 1 

Thus,  we  refer  to  the  ratio  on  the  right-hand  side  of  eqn.  (5)  as  the  fractional  surface  content 
( FSC )  of  the  low-order  representation  containing  the  first  M  modes  as  it  provides  a  measure 
of  the  detail  from  the  full  surface  retained  in  the  A/-order  reconstruction. 

Replicas  of  the  full  surface  as  well  as  the  low-order  models  were  fabricated  using  a  rapid¬ 
prototyping  method  based  on  powder  deposition  with  a  spatial  resolution  of  80  /iin  as  described 
in  Wu  and  Christensen  [35].  The  full  surface  was  fabricated,  in  addition  to  two  low-order  models 
constructed  from  the  most  dominant  modes  of  the  aforementioned  decomposition:  M  —  5  and 
M  =  16  (M  —  5  implies  only  the  first  five  basis  functions  are  included  while  M  =  16  indicates 
that  the  first  sixteen  basis  functions  are  included  in  the  reconstruction).  The  low-order  models 
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Surface 

FSC 

k  (mm) 

ktms  (mm) 

Skewness 

Flatness 

Full  surface 

1 

4.25 

1.0 

0.16 

2.27 

Af  =  16  model 

0.95 

4.09 

0.975 

0.20 

2.34 

Af  =  5  model 

0.71 

3.66 

0.84 

0.24 

2.82 

Table  1:  Various  characteristics  of  the  rough  surfaces  under  consideration. 


of  the  full  surface  are  hereafter  referred  to  by  the  highest-order  mode,  A/,  included  in  the  low- 
order  reconstruction.  Figures  1(b)  and  1(c)  present  topographical  maps  of  the  M  =  16  (fcjJ£L16  = 
0.975  mm)  and  M  =  5  (^7=5  =  0.84  mm)  low-order  representations,  respectively,  which  clearly 
illustrate  the  increased  topographical  detail  embodied  in  the  models  with  increasing  modal 
content.  In  particular,  the  M  =  16  model  [figure  1(b)]  bears  a  striking  resemblance  to  the 
full  surface  [figure  1(a)]  which  includes  all  383  modes.  For  reference,  the  FSC  of  the  M  =  16 
reconstruction  is  0.95  (^7=16  =  0.975fcJJ{f),  meaning  nearly  all  of  the  topographical  detail  of 
the  full  surface  is  embodied  in  this  model.  In  contrast,  the  low-order  representation  containing 
the  first  five  modes  only  loosely  reflects  the  largest-scale  surface  defects  present  in  the  full 
surface  which  is  justified  by  its  relatively  low  FSC  of  0.71  =  0.81  k}™)-  The  average 

characteristics  of  each  roughness  case  are  summarized  in  table  1 . 

Since  the  low-order  surface  representations  presented  in  figures  1(b)  and  1(c)  were  generated 
by  truncating  the  series  representation  in  eqn  (1)  at  mode  Af,  there  exists  a  set  of  residual 
topographical  modes  from  A/  +  1  to  L  that  represent  the  surface  content  excluded  in  a  low- 
order  representation  given  by  Af.  One  can  reconstruct  the  residual  surface  features  embodied 
in  these  excluded  modes  for  a  given  low-order  representation  Af  by  computing  its  associated 
residual  topographical  field,  as 


L 

z )  =  ??fuu(x,  z)  -  7^fw(x,  z)  =  ^2  aj4>j(x ,  z).  (6) 

j  = A/4*  1 


Figures  1(e)  and  1(f)  present  topographical  maps  of  these  residual  fields  for  the  Af  =  16  and 
Af  =  5  low-order  reconstructions.  The  contour  levels  in  these  maps  are  consistent  with  those 
employed  in  the  maps  of  the  low-order  representations  themselves  in  order  to  facilitate  direct 
comparison.  As  one  would  expect,  a  significant  amount  of  the  full-surface  content  is  apparent 
in  the  residual  field  for  the  Af  =  5  case  [figure  1(f)],  as  its  spatial  characteristics  are  quite 
consistent  with  the  overall  character  of  the  full  surface.  This  observation  is  consistent  with  the 
fact  that  the  Af  =  5  low-order  model  contains  only  71%  of  the  full  surface  content,  meaning  that 
the  residual  field  contains  a  not  insignificant  amount  of  information  of  the  full  surface  (29%).  In 
contrast,  inclusion  of  11  additional  modes  in  the  Af  =  16  case  yields  a  residual  topographical  field 
[figure  1(e)]  that  bears  little  resemblance  to  the  full  surface.  Instead,  the  Af  =  16  residual  field 
is  marked  by  weak,  smaller-spatial-scale  fluctuations  in  surface  topography.  The  dominance  of 
smaller-spatial-scale  topographical  fluctuations  in  this  residual  field  is  consistent  with  the  fact 
the  higher-order  modes  generated  by  SVD  embody  the  contributions  of  sinaller-scale  spatial 
features. 

Finally,  figure  1(d)  presents  streamwise  profiles  of  the  full-surface  roughness  as  well  as  the 
two  low-order  representations  coincident  with  the  lines  in  figures  l(a)-(c).  These  lines  de¬ 
marcate  the  spanwisc  position  of  the  streamwise- wall-normal  measurement  plane  for  the  PIV 
experiments  presented  herein.  These  profiles  provide  a  more  quantitative  comparison  of  the 
low-order  surface  models  to  the  full  surface.  This  direct  comparison  highlights  the  finer-scale 
details  that  are  excluded  in  the  low-order  representations,  particularly  in  the  case  of  the  5-mode 
model  which  only  captures  the  larger- wavelength  features  of  the  full  surface.  In  contrast,  the 
16- mode  model  provides  a  much  more  accurate  description  of  the  full  surface,  though  loss  of 
smaller- wavelength  information  in  this  model  is  still  quite  apparent.  In  particular,  the  large- 
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Figure  2:  Probability  density  functions  (pdfs)  of  (a)  roughness  amplitude,  ?/,  about  the  mean  elevation, 
(b)  residual  roughness  amplitude  about  the  mean  elevation,  (c)  streamwise  surface  gradient,  di]/dx<  and 
(d)  spanwise  surface  gradient,  dr)/dz. 


scale  protrusions  of  the  16-mode  model  appear  much  smoother  than  those  of  the  full  surface 
which  tend  to  contain  more  jagged  smaller- scale  detail.  These  similarities  and  differences  are 
also  reflected  in  figure  2(a)  which  presents  probability  density  functions  (pdfs)  of  roughness 
amplitude  about  the  mean  elevation  for  each  surface  ease.  The  low- pass-filtering  effect  of  the 
low-order  reconstructions  is  apparent  when  comparing  the  pdf  of  roughness  amplitude  for  the 
M  =  5  model  with  that  of  the  full  surface.  Removal  of  the  finer  scales  embodied  in  the  higher- 
order  SVD  modes  effectively  smooths  the  full-surface  topography,  though  the  M  =  16  model 
captures  the  overall  topographical  trends  of  the  full  surface  quite  well.  The  pdfs  of  the  residual 
modes  for  the  M  =  5  and  M  =  16  cases  are  presented  in  figure  2(b)  and  confirm  that  significant 
intermediate-  and  fine-scale  detail  of  the  surface  is  lost  in  the  M  =  5  case  (~  ±1.5  mm)  while  the 
pdf  for  the  residual  roughness  amplitudes  for  the  M  =  16  case  is  much  narrower  (~  ±0.5  mm). 
This  filtering  effect  of  truncation  is  most  apparent,  however,  in  pdfs  of  the  streamwise  and 
spanwise  surface  gradients,  drj/dx  and  drj/dz ,  presented  in  figures  2(c)  and  2(d),  respectively. 
Truncation  of  the  higher-order  modes  clearly  reduces  the  spatial  gradients  of  the  topography, 
most  notably  in  the  case  of  the  M  =  5  model,  though  the  M  =  16  model  also  deviates  from  the 
full-surface  behavior  much  more  than  is  discernable  from  the  roughness  amplitude  pdfs.  These 
trends  therefore  highlight  the  fact  that  the  higher-order  SVD  modes  not  only  embody  smaller- 
scale  elements  of  the  full  surface  but  also  contribute  to  the  larger-scale  protrusions  themselves, 
particularly  the  jagged  edges  of  these  protrusions  as  is  notable  in  figure  1(d). 

To  facilitate  testing,  the  aforementioned  printer  was  used  to  construct  replicas  of  each  to¬ 
pography  layer-by-layer  with  a  maximum  spatial  footprint  of  25  x  30  cm2  and  a  mean  thickness 
of  approximately  6  mm.  In  order  to  accommodate  these  roughness  panels,  the  upstream  half 
of  the  suspended  flat  plate  in  the  wind  tunnel  was  raised  relative  to  the  downstream  half  such 
that  the  mean  elevation  of  the  roughness  was  coincident  with  the  upstream  smooth  wrall.  Since 
the  original  spatial  footprint  of  the  digitized  topography  was  not  sufficient  to  fill  this  large  of  an 
area,  the  topography  was  mirrored  in  both  the  streamwise  and  spanwise  directions  to  achieve  an 
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appropriate  streamwise  fetch  of  roughness  [35].  Two  different  streamwise  fetches  of  roughness 
were  considered  in  this  effort:  lm  (equivalently  ~  105)  and  3m  (~  305)  to  facilitate  the  study 
of  both  developing  and  developed  turbulent  boundary  layers.  The  boundary  layers  were  tripped 
to  initiate  transition  to  turbulence  at  the  leading  edge  of  the  suspended  plate  and  grew  over 
smooth-wall  conditions  for  the  first  3  m  and  5  m  for  the  developed  and  developing  flow  exper¬ 
iments,  respectively,  after  which  they  encountered  the  roughness  conditions  considered  herein. 
For  reference,  a  2  m-long  smooth  cast  aluminum  plate  was  laid  along  the  downstream  half  of 
the  boundary-layer  plate  to  extend  the  upstream  smooth-wrall  fetch  an  additional  2  m  for  the 
developing-flow  cases,  after  which  the  flow  encountered  the  1  m-long  patch  of  roughness.  Since 
the  initial  boundary-layer  development  occurred  over  a  smooth  wall,  this  1  m-long  roughness 
scenario  therefore  represents  the  impact  of  short  streamwise  fetches  of  roughness  on  an  incoming 
smooth-wall  turbulent  boundary  layer.  As  such,  one  would  expect  the  formation  of  an  inter¬ 
nal  layer  within  the  turbulent  boundary  layer  wherein  roughness  effects  are  confined.  In  the 
3  m-long  roughness  scenario,  the  flow  is  termed  developed  to  reflect  the  fact  that  the  internal 
roughness  layer  has  grown  to  occupy  the  entire  wall-normal  extent  of  the  turbulent  boundary 
layer  and  a  self-similar  state  is  attained.  The  reproduction  of  over  thirty  individual  roughness 
panels  for  each  case  was  required  to  achieve  a  streamwise  fetch  of  1  m  while  an  additional  sixty 
panels  per  case  were  required  to  extend  the  streamwise  fetch  to  3  m.  The  panels  were  mounted 
to  cast  aluminum  plates  which  were  then  laid  along  the  downstream  half  of  the  wind  tunnel. 
All  measurements  were  made  0.8  m  (~  85)  and  2.5  m  (~  255)  downstream  of  the  leading  edge 
of  the  roughness  for  the  developing-  and  developed-flow  cases,  respectively. 

2.2  Streamwise-wall-normal  plane  PIV  details 

Over  forty-five  hundred  statistically  independent,  two-dimensional  velocity  (u,  v)  fields  were  ac¬ 
quired  by  PIV  in  the  streamwise-wall-normal  (x,  y)  plane  at  the  spanwise  center  of  the  tunnel 
for  each  surface  condition  to  minimize  sampling  errors  in  the  computed  statistics.  It  should 
be  noted  that  the  streamwise  and  spanwise  positions  of  this  measurement  plane  were  carefully 
maintained  over  the  same  roughness  features  for  the  developing-  and  developed-flow  experi¬ 
ments,  respectively,  allowing  a  meaningful  comparison  of  the  flow  statistics.  The  flow  field  was 
illuminated  with  a  500  /zm-thick  laser  sheet  generated  by  a  pair  of  Nd:YAG  lasers  (200  mJ/pulse, 
5  ns  pulse  duration)  and  a  combination  of  spherical  and  cylindrical  lenses.  A  high-energy  mirror 
directed  the  laser  sheet  into  the  wind  tunnel  such  that  it  was  normal  to  the  flow  boundary 
and  parallel  to  the  flow  direction.  The  flow  was  seeded  with  lfim  olive  oil  droplets  and  time- 
separated  images  of  the  scattered  light  from  the  particles  were  captured  with  a  4k  x  2.8k,  12-bit 
frame-straddle  CCD  camera  over  a  field  of  view  of  1.155  x  0.85  (streamwise  by  wall-normal) 
and  1.45  x  5  for  the  developing  and  developed  flow  experiments,  respectively.  The  roughness 
at  the  measurement  location  was  painted  black  to  reduce  reflections  of  laser  light;  however,  the 
remaining  unsuppressed  reflections  rendered  measurements  in  the  region  y  <  0.06  —  0.15  impos¬ 
sible  for  the  rough-wall  cases.  For  reference,  all  rough-wrall  data  presented  herein  is  compared 
to  data  from  a  smooth- wall  PIV  experiment  conducted  at  Re#  =  11400  (5+  =  3350). 

The  pairs  of  PIV  images  were  interrogated  using  a  recursive,  two-frame  cross-correlation 
method.  The  sizes  of  the  interrogation  windows  were  chosen  to  maintain  a  consistent  vector  grid 
spacing  between  the  various  cases  when  scaled  in  inner  units  (i.e.  by  ur  and  v)  and  the  second 
window  was  offset  by  the  bulk  displacement  to  reduce  errors  associated  with  loss  of  image  pairs. 
The  grid  spacing  for  all  cases  is  Ax4-  =  A y^~  =  A4-  ^18  (where  («)+  denotes  normalization 
in  inner  units)  using  the  developed-flow  viscous  length  scales  for  normalization.  The  resulting 
velocity  vector  fields  were  then  validated  using  neighborhood  mean  and  median  comparisons  to 
remove  spurious  vectors.  On  average,  a  valid  vector  yield  of  >  97%  was  achieved,  minimizing  the 
need  for  interpolation  of  holes.  Finally,  each  velocity  field  was  low-pass  filtered  with  a  narrow 
Gaussian  filter  to  remove  noise  associated  with  frequencies  larger  than  the  sampling  frequency 
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Figure  3:  Schematics  for  (a)  developing  and  (b)  developed  flow  PIV  measurements  in  the 
streamwise-wall-normal  plane. 


of  the  interrogation. 

With  regard  to  PIV  measurement  uncertainties,  errors  can  arise  from  multiple  sources. 
In  particular,  the  random  error  associated  with  determining  particle  displacements  in  PIV  is 
approximately  5%  of  the  particle-image  diameter  [37].  In  the  present  study,  the  mean  particle- 
image  diameter  was  approximately  2  pixels,  yielding  a  random  error  of  0.1  pixels.  Therefore, 
since  the  bulk  displacements  in  the  present  experiments  were  10-12  pixels,  this  random  error 
is  less  than  1%  of  the  full-scale  velocity.  In  addition,  bias  errors  can  play  an  important  role 
in  the  accuracy  of  PIV  data.  Bias  due  to  loss  of  image  pairs  was  minimized  by  utilizing  a 
larger  second  interrogation  window  and  a  bulk  window  offset  during  interrogation  of  the  PIV 
images.  Bias  errors  due  to  the  peak-locking  effect  were  also  minimized  as  the  particle-iinagc 
diameters  exceeded  2  pixels  [38,  39],  Finally,  sampling  errors  in  the  statistics  presented  arc 
small  (approximately  1%)  since  averaging  was  performed  over  at  least  1,250,000  statistically 
independent  velocity  samples  per  data  point.  As  such,  symbol  size  throughout  accommodates 
the  uncertainty  bounds  for  each  statistic  presented. 

Measurements  were  conducted  under  both  developing  and  developed  flow  conditions  for  flow 
over  the  original  roughness  and  the  two  low-order  representations.  Figure  3  presents  a  schemat  ic 
of  both  surface  scenarios  whereby  the  developing-flow  measurements  wrere  conducted  0.8  in  (~ 
8 5)  downstream  of  the  abrupt  transition  from  smooth  to  rough  conditions  while  the  developed 
flow  measurements  were  conducted  2.8  m  28(5)  downstream  of  this  surface  transition  where 
self-similar  conditions  are  assured.  In  both  scenarios,  the  flow  is  allowed  to  develop  into  a 
smooth- wall  self-similar  turbulent  boundary  layer  prior  to  reaching  the  rough  surface. 

2.3  Streamwise-spanwise  plane  stereo  PIV  details 

For  the  stereo  PIV  measurements  in  the  x  —  z  plane  at  y  =  0.047(5  (relative  to  the  virtual 
origin  for  each  rough- wall  case),  the  two  cameras  were  positioned  approximately  1.5  m  away 
from  the  measurement  plane  and  imaged  the  flow  through  a  transparent  section  in  the  wind- 
tunnel  ceiling  (see  figure  4).  The  imaging  paths  of  the  cameras  were  rotated  at  ±13°  with 
respect  to  the  wall-normal  ( y )  axis.  The  laser  light  sheets  were  formed  in  a  manner  identical  to 
that  described  earlier  but  were  introduced  into  the  tunnel  through  a  glass  sidewall  while  their 
orientations  were  carefully  adjusted  to  ensure  they  remained  parallel  to  the  wall  at  y  =  0.047(5. 
Uniform  image  focus  was  ensured  in  both  cameras  across  the  entire  field  of  view  by  satisfying 
the  Scheimpflug  condition.  The  pairs  of  images  acquired  by  each  camera  were  interrogated  and 
validated  independently  in  accordance  with  the  methodology  described  above,  resulting  in  pairs 
of  instantaneous  planar  fields  of  two-dimensional  particle  displacements  from  the  two  cameras. 
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Figure  4:  Schematic  for  developed  flow  stereo  PIV  measurements  in  the  streamwise-spanwise 
plane  at  y  —  0.047£. 


Figure  5:  Influence  of  background  reflections  on  single-point  statistics,  (a)  raw  image,  (b)  wall-normal 
Reynolds  normal  stress  <  v'2  >+. 


These  pairs  of  displacement  fields  were  then  reeombined  into  instantaneous  veloeity  fields  using  a 
mapping  function  generated  via  calibration  of  the  imaging  system.  This  calibration  employed  a 
dual-plane  (1-mm  separation)  target  containing  white  dots  spaced  at  in-plane  intervals  of  10  mm 
over  a  20  cm  x  20  em  area.  The  use  of  this  dual-plane  target  alleviated  the  need  to  physieally 
translate  the  target  in  the  depth  direction  through  the  lightsheet  thickness.  Instead,  the  target 
was  earefullv  aligned  in  the  field  of  view  so  that  its  mid-plane  coincided  with  the  laser  lightsheet. 
With  the  target  fixed  at  this  location,  an  image  was  acquired  from  each  camera,  and  these  images 
were  used  to  generate  the  mapping  function  via  the  least-squares  method  of  Soloff  et  ah  [40] 
using  third-order  polynomials  for  the  in-plane  coordinates  and  a  first-order  polynomial  for  the 
out-of-plane  coordinate.  A  draw-baek  of  this  technique  is  that  obtaining  perfeet  alignment 
of  the  target  with  the  laser  lightsheet  is  virtually  impossible.  To  overcome  this  problem,  the 
self-ealibration  seheme  originally  proposed  by  Wieneke  [41]  was  used  to  optimize  the  mapping 
function.  The  final  mapping  function  was  then  used  to  reconstruct  three-dimensional  veloeity 
vectors  on  the  measurement  plane  from  the  pairs  of  two-dimensional  particle  displacements. 

As  these  rough-wrall  experiments  were  begun,  it  was  realized  that  an  important  challenge 
must  be  overcome.  In  particular,  as  the  laser  lightsheet  was  parallel  to  and  wTithin  a  few  millime¬ 
ters  of  the  rough  surfaces,  surface  reflections  beeame  a  debilitating  issue.  While  smooth-wall 
measurements  were  not  tainted  by  sueh  effeets,  the  eomplex  nature  of  the  surfaee  roughness 
under  consideration  led  to  intense  reflections  of  laser  light  that  was  directed  in  all  directions 
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Figure  6:  Image  after  fixed  objects  subtraction  using  original  algorithm  of  Deen  et  al.  [42].  (a)  General 
view;  (b)  Zoom  in  on  a  former  bright  spot. 


and  thus  appeared  in  the  PIV  images  acquired.  To  highlight  this  issue,  a  sample  PIV  image 
for  flow  over  the  M  =  5  model  at  y  =  0.047^  is  shown  in  figure  5  and  the  strong  background 
reflections  due  to  the  largest  roughness  elements  are  obvious,  particularly  how  intense  these  un¬ 
wanted  reflections  are  compared  to  the  imaged  light  from  the  tracer  particles.  These  background 
reflections  unfortunately  corrupt  the  PIV  interrogation  and  induce  strong  bias  in  single-point 
statistics  [(t’'2)4-  shown  in  figure  5(b)]. 

A  solution  to  this  type  of  background  reflection  problem  was  recently  proposed  by  Deen  et 
al.  [42].  The  basic  idea  of  this  method  is  to  subtract  frame  D  from  frame  A  of  instantaneous 
PIV  realizations,  which  in  theory  would  render  the  intensity  of  the  fixed  objects  zero  since  they 
would  scatter  similar  light  in  both  frames  (compared  to  the  random  nature  of  the  scattered  light 
from  the  moving  tracer  particles).  To  compensate  for  intensity  differences  due  to  disparity  in 
the  power  of  the  two  laser  heads,  the  authors  introduced  an  intensity  normalization  based  on 
local  values  of  minimum  and  maximum  intensity.  To  sample  for  extreme  values  and  normalize 
locally,  they  used  an  image  spot  of  size  larger  that  the  average  particle  image  but  smaller  than 
the  intended  PIV  window  size.  The  local  normalization  is  then  carried  out  as 

A.M=  /(X)  /  (X)  (7) 

1minv *) 

where  /(x)  is  the  sliding  local  intensity  and  Imax(x)  and  /min(x)  are  its  maximum  and  minimum 
values,  respectively.  After  normalizing,  since  the  intensity  values  are  bounded  between  0  and  1, 
the  two  PIV  frames  can  be  subtracted  from  each  other.  After  subtracting  frame  A  from  frame 
B ,  frame  A  will  be  recovered  by  extracting  the  positive  values  of  intensity  from  the  result  while 
frame  B  will  be  recovered  by  extracting  the  negative  values  of  intensity  from  the  result  and 
taking  their  absolute  value.  This  approach  worked  well  for  the  application  presented  by  Deen 
et  al.  [42];  however,  as  can  be  observed  in  figure  6b,  the  normalization  spot  leaves  a  footprint  in 
the  formerly  bright  regions  of  the  PIV  images.  Since  the  local  maximum  intensity  may  change 
from  one  frame  to  another  due  to  particle  displacement,  the  relative  intensity  of  the  background 
would  also  be  different  after  local  normalization.  Thus,  as  proposed,  this  correction  method  is 
quite  sensitive  to  extreme  values. 

To  overcome  this  problem  of  extreme  value  sensitivity,  we  modified  the  method  using  an 
initial  local  normalization  based  on  the  local  median  and  minimum  intensity  values  as 

JV(x)  = — _ _  (8) 

[>  /median(x)-/mm(x)-  W 

A  frame  subtraction  was  then  applied  followed  by  a  local  normalization  according  to  equa¬ 
tion  (7).  Finally,  the  intensity  values  were  stretched  to  a  16- bit  intensity  range.  The  final  result 
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Figure  7:  Image  after  modified  fixed  objects  subtraction,  (a)  General  view;  (b)  Zoom  in  on  a  former 
bright  spot. 


of  this  modified  method  is  presented  in  figure  7  wherein  the  effect  of  the  normalization  spot  is 
virtually  non-existent  which  greatly  reduced  any  bias  due  to  background  reflections.  This  mod¬ 
ified  background  subtraction  technique  was  employed  for  all  of  the  stereo  PIV  measurements 
described  herein,  including  the  smooth-wall  PIV  images  for  consistency. 


3  Analysis 


The  efficacy  of  the  two  different  low-order  representations  in  reproducing  the  characteristics 
of  flow  over  the  full  surface  under  both  developing-  and  developed-flow  scenarios  is  evaluated 
by  directly  comparing  various  single-point  statistics,  including  the  mean  velocity,  the  Reynolds 
normal  stresses  ((t/2)  and  (r,/2))  and  the  Reynolds  shear  stress  (RSS;  ( u'v ')).  In  all  cases,  the 
mean  velocity  profile  was  formed  by  ensemble-averaging  all  velocity  realizations  for  a  given  case 
followed  by  line- averaging  in  the  streamwise  direction.  The  turbulence  statistics  considered 
herein  were  then  computed  in  a  similar  manner  using  fluctuating  velocity  fields  derived  from 
the  aforementioned  mean  velocity  profiles. 

Probability  density  functions  are  also  used  to  contrast  the  instantaneous  contributions  to 
the  RSS  for  flow  over  the  two  low-order  models  and  the  full  surface.  In  this  regard,  the  in¬ 
stantaneous  u'v'  events  that  are  averaged  to  generate  the  mean  RSS  profile  can  be  formed  by 
different  combinations  of  v!  and  v'  depending  upon  which  quadrant  of  the  u'  —  v'  plane  a  given 
instantaneous  RSS  event  resides.  In  particular,  negative  contributions  to  (u'v')  are  attributable 
to  ejection  ( Q2 :  v!  <  0,  vf  >  0)  and  sweep  ( Q\ :  v!  >  0,  vf  <  0)  events  while  positive  contributions 
are  generated  by  inward  (Q3:  v!  <  0,  v'  <  0)  and  outward  (Qj:  u'  >  0,  vf  >  0)  interactions. 
To  further  explore  the  efficacy  of  the  low-order  models  in  reproducing  the  character  of  these 
individual  quadrant  events  for  flow  over  the  full  surface,  quadrant  analysis,  as  first  proposed 
by  Lu  and  Willmarth  [43],  is  also  applied  to  all  cases.  In  quadrant  analysis,  the  mean  RSS 
at  each  wall-normal  position  is  decomposed  into  contributions  from  four  quadrants  excluding  a 
hyperbolic  hole  of  size  H  as 

!  p 

{u'v')q(y,  H)  =  -Y^u'{xj,y)v'{xj,y)lQ{xj,y\  H),  (9) 

J  =  1 


where  P  is  the  total  number  of  velocity  vectors  at  each  wall-normal  position  and  Iq  is  the 
indicator  function  defined  as 


iQ(xj,y\ 


1,  when  \u'(xj,y)v'(xj,y)\Q  >  T 
0,  otherwise, 


(10) 
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Surface 

Reg 

Ue  (m/s) 

S  (mm) 

Full  surface 

15970 

16.9 

125.8 

M  =  16  model 

15540 

16.8 

120.9 

M  =  5  model 

15400 

16.9 

126.2 

Table  2:  Relevant  experimental  parameters  for  developing-flow  experiments. 


where  T  is  a  threshold  that  allows  one  to  consider  various  magnitudes  of  instantaneous  RSS 
events  that  contribute  to  the  mean  RSS.  For  developing  flow  this  threshold  is  taken  to  be 
T  =  i7|{uV)|£Jux  (the  maximum  in  the  mean  RSS  for  the  full-surface  case)  while  for  developed 
flow  the  threshold  is  defined  as  T  =  Hcru(y)crv(y)  where  au  =  (u'2)1/2  and  av  =  (r'2)1/2  arc  RMS 
streamwise  and  wall-normal  velocities,  respectively.  Here,  the  value  H  represents  a  threshold 
on  the  strength  of  the  RSS-producing  events  considered  in  the  analysis,  with  H  =  0  allowing 
inclusion  of  all  u!v*  events  and  increasing  values  of  H  allowing  inclusion  of  only  increasingly 
intense  RSS-producing  events.  Using  this  thresholding  of  uV  events,  one  can  also  document 
the  fraction  of  space,  Nq ,  they  occupy  as 

nq{v,h)  =  (ii) 


Finally,  two- point  velocity  correlation  coefficients  are  compared  under  developed  flow  con¬ 
ditions  to  discern  the  impact  of  the  low-order  surface  models  on  the  spatial  structure  of  the 
flow  compared  to  flow  over  the  full  surface  as  well  as  the  smooth- wall  baseline.  In  this  regard, 
previous  studies  of  wall  turbulence  have  established  a  strong  linkage  between  the  average  spa¬ 
tial  characteristics  of  these  correlations  and  the  dominant  spatial  structure  of  the  flow.  For  the 
present  streamwise-wall-normal  PIV  measurements,  two-point  velocity  correlation  coefficients 
of  the  form 


pij(Ax,  y;  yref) 


{u'i{x,yttj)u'j{x  +  Ax,  7/)) 
<7i(yref)<Tj(y) 


(12) 


where  Oi  and  o 3  are  the  root-mean-square  (RMS)  velocities  of  the  ith  and  jth  velocity  com- 
ponents,  are  computed  for  the  streamwise  and  wall-normal  velocity  components  yielding  three 
correlation  coefficients  puu ,  pvv  and  puv.  For  reference,  previous  studies  indicate  fair  similarity 
in  these  correlation  coefficients  for  smooth-  and  rough- wall  flow  outside  the  roughness  sublayer, 
though  a  shortening  in  the  streamwise  extent  of  puu  has  been  reported  for  both  idealized  [44,  21] 
and  irregular  [26]  roughness  that  can  extend  well  beyond  the  roughness  sublayer.  Nonetheless, 
such  similarity  is  consistent  with  Townsend’s  hypothesis  and  indicates  relatively  little  impact 
of  roughness  on  the  spatial  structure  of  the  outer  layer. 


4  Streamwise- Wall-Normal  Plane  Measurements 

4.1  Developing  Flow 

This  section  presents  analysis  of  the  PIV  datasets  acquired  0.8  m  (~  8(5)  downstream  of  the  lead¬ 
ing  edge  of  the  roughness,  meaning  that  roughness  effects  arc  confined  within  an  internal  layer 
initiated  at  the  abrupt  transition  from  smooth-  to  rough-wall  conditions.  Table  2  summarizes 
the  relevant  details  of  these  experiments. 

4.1.1  Mean  velocity 

Figure  8  presents  the  mean  streamwise  velocity  profiles  for  the  various  surfaces  under  developing 
flow  conditions.  A  deficit  in  the  mean  velocity  is  evident  for  the  rough-wall  cases  compared  to 
the  smooth-wall  profile.  The  wall-normal  extent  of  this  deficit,  which  can  be  interpreted  as  an 
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Figure  8:  Mean  velocity  profiles  in  outer  units  (normalized  with  Ue  and  £)  for  developing  flow.  Not  all 
data  points  shown  for  clarity.  Q:  Smooth;  □:  M  =  5  model;  A:  M  =  16  model;  0:  Full  surface. 


Figure  9:  Profiles  of  (a)  Reynolds  normal  stresses,  (u'2)  and  (t/2),  and  (b)  Reynolds  shear  stress,  —  (ttV), 
for  developing  flow.  Not  all  data  points  shown  for  clarity.  0:  Smooth;  □;  M  =  5  model*  A:  M  =  16 
model;  0:  Full  surface. 


internal  layer  within  which  roughness  directly  impacts  the  flow,  is  surface-dependent,  with  the 
full  surface  showing  the  largest  wall-normal  extent  of  velocity  deficit  ( y  <  0.6(5).  This  deficit 
increases  with  increasing  modal  content  in  the  low-order  representations  and  is  largest  for  flow 
over  the  full  surface.  For  example,  the  profile  for  the  M  =  5  model  deviates  least  from  the 
smooth- wall  profile,  indicating  a  reduced  impact  of  roughness  on  the  mean  flow.  In  contrast, 
the  mean  profile  for  the  M  =  16  representation  collapses  quite  well  with  the  profile  for  the 
full  surface,  indicating  that  the  16  most  dominant  topographical  modes  (4.2%  of  the  383  total 
modes  but  95%  of  the  FSC)  are  sufficient  for  capturing  the  impact  of  the  full  surface  on  the 
mean  velocity  for  developing  flow. 

4.1.2  Reynolds  stresses 

Figure  9(a)  presents  profiles  of  the  streamwise  and  wall-normal  Reynolds  normal  stresses,  (u'2) 
and  { v /2),  normalized  by  U 2  for  developing  flow.  Enhancement  of  the  Reynolds  normal  stresses 
is  evident  with  increasing  modal  content,  with  both  ( u' 2)  and  (v'2)  enhanced  by  roughly  40%  in 
the  full-surface  case  relative  to  smooth-wall  flow.  With  regard  to  the  low-order  representations. 
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Figure  10:  Probability  density  functions  (pdfs)  of  instantaneous  u'v'  events  at  (a)  y  =  0.1(5,  (b)  y  = 
0  15(5,  (c)  y  =  0  2(5  and  (d)  y  =  0.25(5  for  developing  flow.  Not  all  data  points  shown  for  clarity. 
O’  Smooth;  □:  M  =  5  model;  A:  M  =  16  model;  0-  Full  surface. 


the  M  =  5  model  shows  only  moderate  enhancement  of  ( v'2 )  and  (u'2)  compared  to  the  smooth- 
wall  baseline.  Consistent  with  the  trends  in  the  mean  velocity,  the  M  =  16  result  shows  strong 
consistency  with  the  full-surface  results  except  near  y  =  0.1(5  where  this  low-order  representation 
produces  slightly  lower  values  of  {ua)  compared  to  the  full  surface.  This  observation  highlights 
the  importance  of  even  finer-scale  topographical  details  in  the  local  flow  behavior  in  the  near¬ 
wall  region.  The  wall-normal  extent  of  the  internal  layer  formed  by  the  roughness  grows  with 
increasing  modal  content,  with  the  internal  layer  for  flow  over  the  full  surface  and  the  M  =  16 
model  extending  to  roughly  y  ^  0.55<5  for  both  {ua)  and  (v*2)  but  to  only  y  ^  0.4(5  for  the 
M  =  5  model.  Beyond  these  wall-normal  locations,  the  rough-wall  profiles  collapse  well  with 
smooth- wall  flow. 

Interestingly,  the  M  —  16  model  reproduces  the  full-surface  mean  RSS,  —  (u'v'),  profile 
[figure  9(b)]  extremely  well  over  the  entire  wall-normal  range  over  which  measurements  were 
possible,  including  replication  of  the  wall-normal  extent  of  the  internal  layer  for  the  full  surface 
(y  ±2  0.55).  This  measure  of  the  internal-layer  thickness  for  both  the  M  =  16  model  and  the  full 
surface  is  slightly  smaller  than  that  noted  from  the  mean  velocity  and  Reynolds  normal  stress 
profiles.  In  contrast,  the  M  =  5  model  yields  much  lower  values  of  mean  RSS  as  well  as  a  thinner 
internal  layer  (y  ic  0.355),  though  the  magnitude  of  the  M  =  5  mean  RSS  is  still  larger  than 


16 


0  0.2  0.4  0.6  0  8  1  0  0  2  0  4  0.6  0  8  1 

y  /  8  >  /  8 

Figure  11:  (a,b)  Quadrant  contributions  to  the  mean  RSS,  {u'z/)q/U^ ,  and  (c,d)  space  fractions,  Nq , 

for  developing  flow  as  a  function  of  wall-normal  position  for  H  =  0.  Not  all  data  points  shown  for  clarity. 

O:  Smooth;  □:  M  =  5  model;  A:  M  =  16  model;  0*  Full  surface. 


0.2 


that  for  the  smooth  wall.  To  further  explore  these  trends,  particularly  the  contributions  of  the 
various  events  that  yield  the  overall  mean  RSS  (ejections,  sweeps,  inward/outward  interactions), 
figure  10  presents  pdfs  of  instantaneous  uV  events  at  various  wall-normal  locations  for  flow  over 
all  surface  cases.  Event  magnitudes  are  normalized  by  the  magnitude  of  the  peak  mean  RSS  for 
the  full-surface  case  in  order  to  directly  compare  the  efficacy  of  the  low-order  surface  models  in 
reproducing  the  flow  behavior  of  the  full  surface.  Consistent  with  the  trends  noted  in  the  mean 
RSS  profiles,  the  M  =  5  model  does  not  reproduce  the  pdf  of  the  full  surface,  most  notably 
close  to  the  wall.  These  differences  are  apparent  in  both  the  negative  and  positive  tails  of  the 
pdfs,  meaning  that  all  four  quadrant  events  are  likely  modified  in  the  M  =  5  model  compared  to 
the  full  surface.  In  contrast,  the  M  =  16  surface  model  produces  pdfs  that  are  quite  consistent 
with  the  character  of  their  full-surface  counterparts.  Far  from  the  wall,  the  pdfs  from  all  three 
surfaces  collapse  toward  the  smooth- wall  trends  as  the  outer-edge  of  the  internal  roughness  layer 
is  approached. 
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4.1.3  Quadrant  analysis 

Extending  the  analysis  afforded  by  pdfs  of  instantaneous  RSS-produeing  events,  quadrant  anal¬ 
ysis  is  employed  to  assess  how  consistent  the  various  quadrant  contributions  to  the  mean  RSS 
profiles  are  amongst  all  surface  conditions.  Figures  11(a)  and  11(b)  present  contributions  of 
the  four  quadrant  events  to  the  mean  RSS  for  a  threshold  of  H  =  0  wherein  all  instantaneous 
u'v'  events  are  included  in  the  decomposition.  Consistent  with  the  mean  RSS  profiles,  ejections 
(Q2)  and  sweeps  ( Q\ )  dominate  over  inward  (Q4)  and  outward  (Qi)  interactions  for  all  eases. 
Of  particular  interest,  while  the  M  =  5  model  fails  to  reproduce  the  quadrant  contributions  of 
the  full-surfaee  flow,  all  four  quadrant  contributions  to  the  mean  RSS  for  the  M  —  16  low-order 
representation  are  quite  consistent  with  those  for  the  full  surface.  Outside  the  internal  layer 
generated  by  roughness  (y  ^  0.4(5  and  0.5(5  for  the  M  =  5  model  and  the  Al  =  16  model/full 
surface,  respectively),  the  quadrant  contributions  eollapse  well  with  the  smooth- wall  results. 
In  contrast  to  the  roughness-indueed  enhancement  noted  in  the  quadrant  contributions  to  the 
mean  RSS,  the  spaee  fractions  of  these  quadrant  events  are  found  to  be  similar  for  all  surface 
eases  [figures  11(e)  and  11(d)].  Thus,  roughness  does  not  appear  to  produee  more  numerous 
RSS-produeing  events  but  instead  tends  to  intensify  the  magnitude  of  these  events  compared 
to  smooth- wall  flow  for  H  =  0. 
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Surface 

Re<? 

Ue  (m/s) 

S  (mm) 

uT  (m/s) 

y*  (H 

<5+ 

AU+ 

K 

ok/S 

Full  surface 

14  340 

17.25 

101.9 

0.785 

19.5 

5  215 

8.2 

218 

107 

0.21 

M  =  16  model 

13  300 

17.2 

98.6 

0.785 

19.8 

4  950 

7.9 

207 

95 

0.21 

M  =  5  model 

12  590 

17.3 

96.6 

0.723 

21.1 

4  580 

5.7 

173 

0.19 

Table  3:  Relevant  experimental  parameters  for  developed-flow  experiments. 


Figures  12(a)  and  12(b)  present  quadrant  contributions  for  the  three  rough  surfaces  and  the 
smooth-wall  baseline  for  a  threshold  defined  by  H  =  4  which  allows  only  the  most  intense  uV 
events  to  be  included  in  the  decomposition.  Note  that  only  ejections  ( Q2 )  and  sweeps  ( Q4 ) 
are  presented  since  the  contributions  of  inward  and  outward  interactions  are  near  zero  for  this 
threshold.  Consistent  with  the  H  =  0  trends,  the  M  =  16  surface  model  reproduces  the  full- 
surface  trends  quite  well  while  the  M  =  5  model  does  not.  Of  interest,  while  the  space  fractions 
for  the  H  =  0  threshold  showed  little  difference  in  the  number  of  quadrant  events  generated  by 
the  two  low-order  models  compared  to  the  full-surface  and  smooth-wall  flow,  the  space  fractions 
for  H  =  4  [figures  12(c)  and  12(d)],  which  includes  only  the  most  intense  RSS-producing  events, 
show  clear  differences  from  smooth-wall  flow.  In  particular,  the  M  =  16  model  reproduces  the 
space  fractions  of  both  ejections  and  sweeps  from  the  full-surface  flow  quite  well.  In  contrast,  the 
M  =  5  model  fails  to  reproduce  the  quantitative  behavior  of  the  full-surface  flow.  Nevertheless, 
compared  to  smooth- wall  flow,  all  of  the  rough- wall  cases  are  found  to  generate  more  numerous 
intense  ejection  and  sweep  events  by  a  factor  of  3-6.  Finally,  outside  the  internal  layer  the 
rough-wall  quadrant  contributions  and  associated  space  fractions  again  collapse  well  with  the 
smooth-wall  results  as  should  be  expected  for  developing  flow. 

4.2  Developed  Flow 

This  section  presents  results  from  the  PIV  measurements  made  ~  25 6  downstream  of  the  leading 
edge  of  the  roughness  such  that  the  rough-wall  flows  have  attained  self-similar  states  based  on 
past  measurements  over  such  surfaces  [26].  Table  3  summarizes  relevant  experiment  parameters 
for  these  developed-flow  measurements. 

4.2.1  Mean  velocity 

Figure  13  presents  mean  velocity  profiles  for  flow  over  the  two  low-order  models,  the  full  surface 
and  the  smooth-wall  baseline  in  both  inner  units  (normalized  by  uT  and  v)  and  in  velocity  defect 
form  (normalized  by  uT  and  <5).  The  friction  velocity,  uT,  for  each  case  was  determined  using 
the  total  shear  stress  method  which  assumes  a  region  of  constant  shear  stress  equal  to  the  wall 
shear  stress  in  the  overlap  and  inner  region  of  the  boundary  layer  [19,  26].  These  values  of  ur 
were  then  used  to  determine  the  virtual  origin,  yQ ,  and  the  roughness  function,  A C/+,  for  the 
rough-wall  cases  by  fitting  the  mean  velocity  profile  to  the  expected  logarithmic  profile  in  the 
log  layer  given  by 

U+  =  -\n(y^-yt)  +  A-AU-,  (13) 

tt 

where  k  —  0.41  and  A  =  5.3  are  the  log-law  constants.  Knowledge  of  A enables  one  to  relate 
the  roughness  studied  herein  to  the  sand-grain  experiments  of  Nikuradse[4]  via  an  equivalent 
sand-grain  height,  that  yields  the  same  A as  the  present  rough  surfaces  through  the 
fully-rough  asymptote  given  by 

AU+  =  -  ln(A^")  +  A  —  8.5.  (14) 

It  should  be  noted  that  the  constant  stress  method  for  determining  uT  has  an  uncertainty  of 
approximately  4-6%  [19,  18]  so  the  symbol  size  of  the  statistics  presented  in  this  section  reflect 
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Figure  13:  Mean  velocity  profiles  in  (a)  inner  units  and  (b)  velocity  defect  form  for  developed  flow. 
Solid,  dashed  and  dotted  lines  in  (b)  demarcate  the  5 k/S  positions  for  the  full,  M  =  16  and  M  =  5  cases, 
respectively,  and  not  all  data  points  shown  for  clarity.  0:  Smooth;  D:  M  =  5  model;  A:  M  =  16  model; 
0*.  Full  surface. 


this  uncertainty  bound. 

The  inner-scaled  profiles  [figure  13(a)]  reveal  the  downward  shift  of  the  rough-wall  profiles 
relative  to  the  smooth-wall  profile  due  to  the  increased  skin  friction  at  the  wall.  Of  interest, 
while  the  inner-scaled  mean  profile  for  the  M  =  5  case  resides  between  the  smooth  and  full- 
surfaee  profiles  (A f/+  =  5.7  for  the  M  =  5  model),  the  M  =  16  ease  (A U+  =  7.9)  sits 
just  above,  but  nearly  reproduces,  the  full-surface  profile  (AJ7+  =  8.2).  One  can  therefore 
infer  that  the  largest  topographical  scales  of  the  full  surface  dominate  A U+  for  the  roughness 
considered  herein,  though  a  slight  difference  does  exist  which  could  represent  contributions  of 
excluded  topographical  modes  to  A U+ .  However,  this  possibility  cannot  be  confirmed  as  the 
difference  between  the  16-mode  and  full-surface  A U+  is  within  the  uncertainty  in  determining 
this  quantity  (3-5%).  In  contrast,  the  mean  velocity  profiles  in  defect  form  [figure  13(b)]  show 
strong  consistency  between  all  surface  cases  outside  the  roughness  sublayer  (y  >  0.2(5),  consistent 
with  Townsend’s  wall  similarity  hypothesis,  and  also  confirm  the  attainment  of  a  self-similar 
state  for  all  cases.  Within  the  roughness  sublayer,  all  three  roughness  cases  sit  below  the 
smooth-wall  case  (consistent  with  many  past  studies  of  both  idealized  [13,  14,  15,  19]  and 
realistic  roughness  [25,  26])  yet  both  low-order  representations  reproduce  the  full-surface  profile 
well. 

4.2.2  Reynolds  stresses 

Figure  14  presents  profiles  of  ( uf2)  + ,  (t/2)4  and  — (wV)+  for  all  surface  cases.  Consistent 
with  Townsend's  wall  similarity  hypothesis,  the  smooth,  full-surfaee  and  model  profiles  collapse 
outside  the  roughness  sublayer  ( y  >  0.2 S  =  5k).  However,  some  differences  are  noted  within 
the  roughness  sublayer.  In  the  case  of  ( u‘ r2)+  [figure  14(a)],  all  three  rough- wall  cases  sit  below 
the  smooth-wall  profile  which  is  a  well-known  effect  of  roughness  in  wall  turbulence.  However, 
while  the  M  =  5  and  M  =  16  cases  collapse  with  one-another,  they  both  yield  slightly  lower 
values  of  (it'2)4'  compared  to  the  full  surface.  This  discrepancy  is  not  entirely  unexpected 
since  the  details  of  the  roughness  can  have  a  significant  impact  on  the  local  flow  behavior  close 
to  the  surface.  As  such,  the  higher-order  topographical  modes  missing  in  the  model  surfaces 
appear  to  play  a  measurable  role  in  (u,2)  +  within  the  roughness  sublayer.  With  regard  to  (t/2)4 
[figure  14(a)],  much  better  consistency  is  noted  as  both  the  5-  and  16- mode  models  reproduce  the 
full-surface  behavior  well,  even  in  the  roughness  sublayer.  In  contrast,  both  low-order  surface 
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Figure  14:  Profiles  of  inner-scaled  (a)  Reynolds  normal  stresses,  (u/2)  and  {v'2),  and  (b)  Reynolds  shear 
stress,  —  (uV),  for  developed  flow.  Lines  as  in  figure  13(b)  and  not  all  data  points  shown  for  clarity. 
O’  Smooth;  □:  M  =  5  model;  A:  M  =  16  model;  0*  Full  surface. 
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models  generate  mean  RSS  profiles  [figure  14(b)]  that  mimic  the  overall  full-surface  behavior 
well  even  within  the  roughness  sublayer. 

While  the  mean  RSS  profiles  for  the  low-order  surface  models  seem  to  reproduce  the  mean 
RSS  behavior  of  flow  over  the  full  surface,  this  collapse  need  not  require  that  the  distributions 
of  the  instantaneous  u'v '  events  contributing  to  these  mean  profiles  be  identical  for  flow  over 
the  full  surface  and  the  two  low-order  models.  To  explore  such  issues,  figure  15  presents  pdfs  of 
instantaneous  u'v'  events  contributing  to  the  mean  RSS  profiles  at  y/S  —  0.1,  0.15,  0.2  and  0.3 
computed  for  all  surface  cases.  While  the  pdfs  of  u'v'  events  away  from  the  wall  [figures  15(d)] 
show  strong  similarity  regardless  of  surface  condition,  the  pdfs  of  u'v '  events  within  the  roughness 
sublayer  [figures  15(a)-(c)]  indicate  that  flow  over  both  the  M  =  5  and  M  =  16  surface  models 
do  not  reproduce  the  number  of  intense  and  negative  u'v'  events  observed  for  flow  over  the 
full  surface.  As  such,  the  low-order  surface  models  appear  to  produce  slightly  fewer  extremely 
intense  ejection  and/or  sweep  events  compared  to  the  full  surface.  In  contrast,  the  positive  tails 
of  the  pdfs  in  figure  15  show  better  collapse,  even  in  the  roughness  sublayer,  indicating  that 
these  positive  contributions  to  ( u'v')+  are  less  affected  by  the  finer-scale  details  of  the  surface 
topography. 

4.2.3  Quadrant  analysis 

To  study  these  trends  in  RSS-producing  events  further,  quadrant  analysis  is  employed  to  dis¬ 
tinguish  ejection  from  sweep  contributions  and  inward-  from  outward-interaction  contributions 
to  the  mean  RSS  profile  for  all  surface  cases.  Figures  16(a)  and  16(b)  present  quadrant  contri¬ 
butions  for  a  threshold  given  by  H  =  0,  meaning  that  all  contributions,  weak  through  intense, 
are  included  in  the  decomposition  As  should  be  expected  from  the  behavior  of  the  mean  RSS 
profiles  [figure  14(b)]  and  the  pdfs  of  u'v'  events  (figure  15),  Q2  and  Q4  contributions  dom¬ 
inate  over  Q\  and  Q3  contributions.  In  particular,  the  profiles  of  Q2  and  Q4  contributions 
are  consistently  3-4  times  larger  in  magnitude  than  those  of  Q\  and  Q3  events.  Further,  the 
rough-wall  results  show  strong  consistency  with  the  smooth-wall  baseline  outside  the  roughness 
sublayer,  again  consistent  with  Townsend’s  hypothesis.  With  regard  to  the  consistency  between 
the  profiles  of  the  low-order  models  and  those  of  the  full  surface,  while  the  Q\,  Q3  and  Q4  pro¬ 
files  collapse  regardless  of  surface,  the  contributions  of  ejections  (Q2)  show  surface  dependence 
close  to  the  wall  with  both  surface  models  yielding  slightly  weaker  contributions  of  ejections  to 
the  mean  RSS  for  the  case  of  H  =  0.  This  behavior  is  consistent  with  the  surface-dependent 
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Figure  15:  Probability  density  functions  (pdfs)  of  instantaneous  inner-scaled  u'v'  events  at  (a)  y  —  0.16, 
(b)  y  =  0.156,  (c)  y  =  0.26  and  (d)  y  =  0.36  for  developed  flow.  Not  all  data  points  shown  for  clarity. 
O’  Smooth;  □:  M  =  5  model;  A:  M  =  16  model;  0:  Pull  surface. 
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behavior  noted  in  the  pdfs  of  ufvf  events  where  the  full  surface  was  found  to  produce  slightly 
more  numerous  intense  negative  uV  events  than  the  two  low-order  surface  models.  As  such, 
this  observation  may  indicate  that  it  is  more  numerous  ejection  events  in  the  full-surface  case 
that  are  primarily  responsible  for  the  surface-dependent  behavior  noted  in  the  pdfs  of  figure  15 
within  the  roughness  sublayer.  Despite  this  difference,  the  fraction  of  space,  Nq,  occupied  by 
these  various  events  [figures  16(a)  and  16(b)]  are  quite  similar  for  the  different  surface  cases 
except  near  the  boundary-layer  edge  where  the  weak  magnitude  of  the  velocity  fluctuations  can 
introduce  increased  uncertainty  in  quadrant  decomposition. 

Figure  17  presents  quadrant  contributions  and  space  fractions  for  a  threshold  of  H  =  4, 
meaning  that  only  the  contributions  of  the  most  intense  tiV  events  to  the  mean  RSS  are 
considered.  While  the  M  =  5  and  M  =  16  surface  models  clearly  yield  weaker  contributions 
of  intense  Q2  events  compared  to  the  full-surface  case  (larger  discrepancy  but  still  qualitatively 
consistent  with  the  H  =  0  trends),  these  surfaces  actually  generate  slightly  higher  contributions 
from  Q4  events  compared  to  the  full-surface  result.  As  such,  the  finer-scale  topographical 
details  of  the  full  surface  excluded  from  the  low-order  surface  models  may  actually  contribute 
in  a  meaningful  way  to  the  generation  of  intense  Q2  and  Q 4  events.  Despite  these  near- wall 
differences,  good  consistency  is  noted  between  the  smooth-  and  rough-wall  results  outside  the 
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Figure  16:  (a,b)  Quadrant  contributions  to  the  mean  RSS,  (uV)q,  and  (c,d)  space  fractions,  Aq,  for 
developed  flow  as  a  function  of  wall-normal  position  for  //  =  0.  Lines  as  in  figure  13(b)  and  not  all  data 
points  shown  for  clarity.  Q:  Smooth;  □:  M  =  5  model;  A:  M  =  16  model;  0:  Full  surface. 


roughness  sublayer  in  accordance  with  outer-layer  similarity. 

4.2.4  Two-point  velocity  correlations 

Figure  18  presents  two-point  correlation  coefficients  of  streamwise  velocity,  puu ,  at  y/6  =  0.1, 
0.15,  0.2  and  0.3  for  flow  over  all  surface  conditions.  While  the  correlations  are  computed 
within  the  two-dimensional  measurement  plane,  they  are  presented  in  one-dimensional  form 
herein  to  quantitatively  assess  their  dependence  on  streamwise  separation,  At.  Well  within 
the  roughness  sublayer  at  y  =  0.15  [figure  18(a)],  while  the  correlations  at  small  At  collapse 
irrespective  of  surface  condition,  there  exists  a  shortening  in  the  streamwise  extent  of  puu  at 
larger  At  for  the  three  rough-wall  cases  is  evident  compared  to  the  smooth-wall  baseline.  These 
trends  indicate  that  the  larger  spatial  scales  of  the  flow  are  more  sensitive  to  roughness  effects 
than  are  the  smaller  scales  of  the  flow.  Interestingly,  the  two  low-order  representations  produce 
a  slightly  more  pronounced  streamwise  shortening  of  puu  than  the  full  surface.  With  increasing 
distance  from  the  wall  this  difference  between  the  full-surface  and  low-order-representation 
profiles  diminishes  quickly  as  the  two  low-order  models  collapse  well  with  the  full-surface  result 
at  y  =  0.155  [figure  18(b)]  and  y  —  0.25  [figure  18(c)],  though  all  three  rough- wall  results 
still  show  shortening  compared  to  the  smooth-wall  baseline.  Outside  the  roughness  sublayer  at 
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Figure  17:  As  in  figure  16  but  //  =  4. 


y  =  0.36  [figure  18(d)],  however,  the  rough- wall  correlations  collapse  well  with  the  smooth-wall 
result  which  is  consistent  with  Townsend's  wall  similarity  hypothesis.  This  collapse  therefore 
suggests  similarity  in  the  average  spatial  structure  of  the  smooth- wall  flow  with  that  of  the 
rough- wall  flows. 

The  two-point  correlation  coefficient  of  wall-normal  velocity,  pvv,  presented  in  figure  19 
shows  less  influence  of  roughness,  though  a  slight  broadening  of  pvv  is  noted  close  to  the  wall  at 
y  =  0.16  [figure  19(a)]  and  y  =  0.156  [figure  19(b)]  for  the  three  rough- wall  cases  compared  to  the 
smooth-wall  baseline.  However,  in  contrast  to  the  trends  noted  at  these  wall-normal  locations 
for  puu,  the  low-order  representations  yield  pvv  profiles  that  match  that  of  the  full-surface  flow 
very  well.  With  increasing  distance  from  the  wall,  pvv  is  found  to  collapse  irrespective  of  surface 
condition.  Comparing  the  streamwise  extent  of  pvv  to  that  of  puu,  it  is  clear  that  the  latter 
is  more  representative  of  the  larger  spatial  scales  of  the  flow  while  the  former  is  more  strongly 
influenced  by  the  smaller  spatial  scales.  Finally,  figure  20  presents  streamwise  profiles  of  the 
cross-correlation  coefficient  of  the  streamwise  and  wall-normal  velocities,  puv.  Similar  to  pvv, 
the  streamwise  extent  of  puv  is  enhanced  within  the  roughness  sublayer  at  y  —  0.16  and  0.156 
(figures  20(a)  and  20(b),  respectively)  compared  to  smooth- wall  flow,  though  the  profiles  for 
the  low-order  representations  collapse  well  with  the  full-surface  results.  With  increasing  wall- 
normal  position  [figures  20(c)  and  20(d)]  these  differences  between  the  smooth-  and  rough-wall 
cases  diminish  indicating  little  impact  of  roughness  on  the  average  spatial  structure  embodied 
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Figure  18:  Streamwise  profiles  of  puu  at  (a)  y  =  0.15,  (b)  y  =  0  155,  (c)  y  =  0  25  and  (d)  y  =  0.35.  Not 
all  data  points  shown  for  clarity.  Q:  Smooth;  □:  M  =  5  model;  A:  M  =  16  model;  0:  Full  surface. 
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in  puv  outside  the  roughness  sublayer. 

5  Streamwise-Spanwise  Plane  Measurements  at  y  —  0.047(5 

Given  the  observations  discussed  above,  particularly  the  notable  differences  observed  between 
flow  over  the  low-order  models  and  the  full  surface  in  the  roughness  sublayer,  wall-parallel 
stereo  PIV  measurements  were  performed  at  y  —  0.0475  for  all  three  rough  surfaces  (y  =  0  in 
the  rough-wall  cases  is  positioned  at  the  virtual  origin  determined  from  the  modified  Clauser 
chart  analysis)  as  well  as  the  smooth- wall  baseline  to  explore  these  differences  further.  Figure  21 
presents  contour  maps  of  the  roughness  topographies  directly  beneath  the  x  —  z  measurement 
planes  for  each  surface  so  that  the  local  topographical  characteristics  can  be  contrasted  with 
the  data  to  be  presented. 

Figures  22  28  present  contour  plots  of  various  single-point  statistics  for  all  four  surface 
cases  with  green  contours  in  the  background  highlighting  the  roughness  patterns  beneath  the 
measurement  plane.  However,  before  discussing  the  details  of  these  fields,  table  4  presents 
the  area- averaged  values  of  the  streamwise  velocity  defect  as  well  as  the  streamwise  and  wall- 
normal  Reynolds  normal  stresses  and  the  RSS.  These  values  are  close,  but  not  identical,  to  their 
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Figure  19:  As  in  figure  18  but  pvv. 


counterparts  from  the  x  —  y  experiments  shown  in  figures  13  and  14.  These  small  differences 
can  be  explained  by  recalling  that  the  x  —  y  measurements  were  performed  at  a  single  spanwise 
location  (demarcated  in  figure  21  by  the  dashed  lines)  while  the  values  presented  in  figure  4  were 
averaged  over  an  area  spanning  many  different  topographical  details  in  the  spanwise  direction. 
Thus,  the  local  details  of  the  topography  can  lead  to  slight  differences  in  the  values  of  the 
single-point  statistics. 

5.1  Mean  velocity 

Prior  to  discussing  the  details  of  the  ensemble-averaged  turbulence  statistics,  it  should  be  noted 
that  the  smooth-wall  results  are  presented  as  a  baseline  for  comparison  as  these  statistics  should 
be  nearly  homogeneous  across  the  (5-scale  field  of  view.  Thus,  any  variability  noted  in  these 
statistics  represents  a  measure  of  their  sampling  error.  Therefore,  any  spatial  variability  ob¬ 
served  in  the  rough-wall  statistics  can  be  compared  to  the  variability  in  the  smooth- wall  baseline 
as  a  means  of  discerning  whether  it  is  due  to  roughness  effects  or  simply  an  artifact  of  sam¬ 
pling.  The  ensemble  average  of  streamwise  velocity  in  defect  scaling,  (Ue  —  u)  +  ,  is  shown  in 
figure  22  for  the  smooth-wall  base  line  as  well  as  the  three  rough-wall  cases  ( M  =  5,  M  —  16 
and  full  surface).  While  the  smooth- wall  result  is  quite  uniform,  all  rough- wall  flows  exhibit 
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Figure  20:  As  in  figure  18  but  puv. 


a  spanwuse-localized  low-momentum  region,  bounded  by  high-momentum  regions.  This  result 
suggests  that  the  roughness  under  consideration  induces  a  “channeling”  effect  in  the  flow.  This 
low-momentum  region  is  located  near  2  ^  0.35  for  the  M  =  5  case  [figure  22(b)],  while  the 
M  =  16  and  full-surface  cases  show  a  low-momentum  region  near  2  ^  0.55  [figures  22(c)  and 
22(d),  respectively].  Further,  it  is  observed  that  the  peak  values  of  the  lowr-  and  high-momentum 
regions  are  higher  in  the  M  —  5  model  than  in  the  other  two  rough-wall  cases. 

The  ensemble-averaged  wall-normal  velocity,  (v)+,  is  shown  in  figure  23.  The  smooth- wall 
result  [figure  23(a)]  is  nearly  zero  across  the  field  of  view,  as  one  would  expect.  However,  while 
the  streamwise  velocity  defect  showed  larger-scale  influences  of  the  roughness,  (v)4  exhibits 
more  localized  influences  of  the  roughness.  In  particular,  (i;)+  is  positive  just  upstream  of  a 
large-scale  peak  in  roughness  and  negative  just  downstream  of  such  a  surface  feature.  Thus,  as 
one  might  expect,  the  flow  is  diverted  away  and  toward  the  wall  by  the  roughness  features  that 
occur  directly  beneath  the  measurement  plane.  These  vertical  excursions  of  the  flow  toward 
and  away  from  the  wall  are  strongest  for  the  M  =  5  case  [figure  23(b)]  which  only  embodies  the 
largest  topographical  scales  of  the  full  surface.  Of  interest,  the  M  =  16  result  [figure  23(c)]  is 
remarkably  consistent  with  the  full-surface  result  [figure  23(d)],  indicating  that  the  intermediate 
topographical  scales  embodied  in  modes  6  16  that  are  missing  in  the  M  =  5  model  play  an 
important  role  in  this  mean  vertical  motion  of  the  flow\  As  mentioned  earlier,  these  intermediate 
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Figure  21:  Local  roughness  topography  below  the  y  =  0.047(5  wall-parallel  fields  of  view  (demarcated 
by  boxed  regions)  for  the  (a)  M  =  5  model,  (b)  M  =  16  model  and  (c)  full  surface.  The  dashed  lines 
represent  the  position  of  the  x  —  y  PIV  measurements,  (d)  Schematic  illustrating  stereo  PIV  view  of  the 
flow  in  the  wall-parallel  field  of  view  above  the  surfaces. 


topographical  scales  add  to  the  jagged  and  irregular  character  of  the  largest  roughness  elements 
present  in  the  topography. 

5.2  Reynolds  stresses 

The  results  for  ensemble- averaged  streamwise  Reynolds  normal  stress,  (t/2)+,  are  presented  in 
figure  24.  As  with  the  mean  velocity  components,  the  smooth-wall  baseline  [figure  24(a)]  shows 
only  slight  variability  across  the  field  of  view'  (~  10%)  which  is  a  measure  of  the  sampling  error 
in  this  second-order  statistic.  The  rough-wrall  results,  however,  show  significant  variability.  Of 
particular  interest,  the  M  =  16  and  full-surface  results  [figures  24(c)  and  24(d),  respectively] 
display  distinct  streamwise-elongated  regions  of  intense  (u/2)+  (one  at  2  «  0.35(5  and  the  other 
at  z  as  0.65(5)  between  which  a  region  of  slightly  weaker  (?i'2)  +  is  noted.  This  elongated  region  of 
weaker  {ua)Jr  is  coincident  with  the  low-momentum  region  noted  in  the  M  =  16  and  full-surface 
streamwise  velocity  defect  results  (figure  22)  while  the  elongated  regions  of  intense  {u'2)+  are 
coincident  with  the  spanwise  boundaries  of  the  low-momentum  region.  A  similar  pattern  is 
also  notable  in  the  M  =  5  {ua)  [figure  24(b)]  where  a  low-momentum  region  was  identified  in 
the  streamwise  velocity  defect  near  2  ~  0.3(5.  In  contrast,  the  ensemble-averaged  wall-normal 
Reynolds  normal  stress,  (u/2)+  (figure  25),  and  the  spanwise  Reynolds  normal  stress,  (u/2)4' 
(figure  26),  show  much  less  large-scale  coherence,  indicating  that  their  behavior  is  more  locally 
governed  by  smaller  spatial  scales  of  the  flow.  Of  interest,  the  M  =  5  model  appears  to  yield 
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Table  4:  Ensemble-  and  area-averaged  values  of  single-point  statistics  at  y  =  0.0475. 


Surface 

(Ue  ~U)  + 

K>+ 

(ttV)+ 

Full 

10.7 

4.6 

1.4 

-0.96 

M  =  16 

10.3 

4.8 

1.4 

-1.00 

M  =  o 

10.0 

4.9 

1.3 

-0.98 

Smooth 

11.8 

4.7 

1.7 

-0.88 

Figure  22:  Contour  maps  of  ensemble- averaged  streamwise  velocity  in  defect  scaling,  ( Ue  —  u)+, 
at  y  =  0.0475.  (a)  Smooth;  (b)  M  =  5  model;  (c)  M  =  16  model;  (d)  Full  surface.  Green  back¬ 
ground  contours  highlight  local  roughness  topography  beneath  the  wall-parallel  measurement 
plane. 


slightly  larger  magnitudes  of  both  (v/2)~h  and  ( wr2)Jr  compared  to  the  M  =  16  model  and  the 
full-surface  results.  In  contrast,  (t,/2)+  and  (it/2)"1"  for  the  M  =  16  model  and  the  full  surface 
match  one-another  well  both  in  magnitude  as  well  as  in  spatial  pattern. 

A  large-scale  spatial  pattern  similar  to  that  observed  in  the  (u/2)+  is  also  noted  in  the 
ensemble- aver  aged  Reynolds  shear  stress,  (uV)+  (figure  27),  where  streamwise-elongated  re¬ 
gions  of  intense  (uV)4*  sit  near  2  «  0.355  and  2  «  0.655  in  the  M  =  16  and  full-surface  results 
while  a  region  of  weaker  (uV)4"  is  observed  in-between.  This  region  of  weaker  (u'vf)+  occurs 
spatially  coincident  to  the  low-momentum  region  noted  in  (Ue  —  u)+  while  the  two  intense  re¬ 
gions  of  (uV)+  lie  along  the  boundaries  of  this  low-momentum  region.  This  imprint  is  also 
apparent  in  the  ensemble- aver  aged  Reynolds  shear  stress  component  (ufwf)+  as  presented  in 
figure  28.  While  the  smooth-wall  baseline  [figure  28(a)]  is  nearly  zero,  as  one  would  expect, 
(ufwf)+  is  found  to  be  intense  in  all  three  rough- wall  cases  with  streamwise-elongated  regions 
of  positive  and  negative  occurring  along  the  boundaries  of  the  low-momentum  regions 

identified  in  (Ue  —  u)  +  .  For  example,  the  M  =  16  and  full-surface  (u/u’/)+  are  positive  in  a 
streamwise-elongated  region  near  z  «  0.355  and  negative  near  z  %  0.65 5.  These  streamwise- 
elongated  regions  occur  spatially  coincident  to  similar  elongated  regions  noted  in  (u'2)'f  and 
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Figure  23:  Contour  maps  of  ensemble-averaged  wall- normal  velocity,  (v)"1",  at  y  =  0.047(5. 
(a)  Smooth;  (b)  M  =  5  model;  (c)  M  =  16  model;  (d)  Full  surface.  Green  background  contours 
highlight  local  roughness  topography  beneath  the  wall-parallel  measurement  plane. 


(uV)+ 

5.3  Quadrant  analysis 

As  discussed  in  section  3  and  as  applied  to  the  PIV  data  in  the  streamwise~wall-normal  plane 
(figures  11,  12,  16  and  17),  quadrant  analysis  can  be  utilized  to  uncover  the  dominant  Revnolds- 
shear-stress- producing  patterns  for  flow  over  the  two  low-order  models  and  the  full  surface  from 
the  streamwise-spanwise  plane  data  deep  within  the  roughness  sublayer.  For  brevity,  only 
the  results  for  a  threshold  of  H  =  4  are  shown  in  figure  29  for  the  M  =  5  and  M  =  16 
models  as  well  as  the  full  surface  for  ejection  (Q2)  and  sweep  (Q4)  events.  Unlike  the  mean 
RSS  in  figure  27  wherein  streamwise-elongated  regions  of  (uV)  are  readily  apparent  along  the 
boundaries  of  the  low-momentum  regions  identified  in  the  mean  velocity  defect  plots  of  figure  22, 
no  such  large-scale  features  are  evident  in  either  the  ejection  or  sweep  contributions  to  the  mean 
Reynolds  shear  stress  for  intense  events  corresponding  to  a  threshold  of  H  =  4.  Thus,  these 
intense  Reynolds-shear-stress-producing  events  are  likely  due  to  snialler-scale  events  that  are 
intimately  tied  to  the  local  roughness  topography. 

5.4  Vortical  activity 

Given  the  streamwise-elongated  imprints  noted  in  the  mean  velocity  defect  {(Ue  —  u)^ ;  figure  22) 
as  well  as  the  streamwise  Reynolds  normal  stress  ((z/2)4-;  figure  24)  and  the  Reynolds  shear 
stress  ({uV)  +  ;  figure  27),  it  is  of  interest  to  assess  whether  enhanced  vortical  activity  might 
be  responsible  for  these  increased  turbulence  levels  in  the  rough-wall  flows.  To  make  such  an 
assessment,  a  vortex  identifier  termed  swirling  strength,  A a-,  is  utilized  to  identify  local  rotational 
motions  in  the  instantaneous  velocity  fields  for  each  surface  condition  in  the  streamwise-  spanwise 
measurement  plane.  Swirling  strength  is  the  imaginary  portion  of  the  complex  eigenvalue  of  the 
local  velocity  gradient  tensor  and  is  an  unambiguous  marker  of  rotation  [45,  46].  Unfortunately, 
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Figure  24:  Contour  maps  of  ensemble- averaged  streamwise  Reynolds  normal  stress,  at 

y  =  0.047£.  (a)  Smooth;  (b)  M  —  5  model;  (c)  M  =  16  model;  (d)  Full  surface.  Green  background 
contours  highlight  local  roughness  topography  beneath  the  wall-parallel  measurement  plane. 


Act  does  not  embody  the  sense  of  the  rotation  as  the  complex  eigenvalues  occur  in  complex- 
conjugate  pairs,  so  A^  is  assigned  the  sign  of  the  local  wall-normal  vorticity  as  XdWy/\wy\  so 
that  clockwise- rotating  vortices  (coy  <  0)  in  the  streamwise-spanwise  measurement  plane  can 
be  distinguished  from  counter-clockwise-rotating  vortices  ( ujy  >  0). 

Figure  30  presents  contour  maps  of  ensemble- averaged  \dujy/\ujy  \  for  the  M  =  5  and  M  =  16 
models  as  well  as  the  full  surface.  While  the  ensemble-averaged  X^y/^y]  for  smooth- wall  flow 
is  essentially  zero  (not  shown  for  brevity),  indicative  of  no  preferential  locations  for  clockwise- 
and  counter-clockwise-rotating  wall-normal  vortices,  it  displays  clear  preferential  locations  for 
such  structures  in  all  three  rough-wall  cases.  Focusing  on  the  M  =  16  model  and  the  full 
surface,  recall  that  a  clear  low-momentum  region  existed  near  2  ^  0.5 S.  The  ensemble- averaged 
An^y/I^yl  for  these  two  cases  display  elongated  regions  of  XriWy/\u)y\  <  0  above  and  XriUy/\u)y\  > 
0  below  this  location  with  a  region  of  \ciUy/\u)y\  «  0  between.  This  pattern  in  both  the  Al  =  16 
model  and  full  surface  is  consistent  with  a  preferential  alignment  of  clockwise-rotating  wall- 
normal  vortices  along  the  upper  boundary  of  the  low-momentum  region  as  well  a  preferential 
alignment  of  counter-clockwise-rotating  wall-normal  vortices  below  this  low-momentum  region. 
The  combined  induction  of  such  vortical  alignments  will  generate  flow  against  the  mean-flow 
direction  and  thus  is  consistent  with  the  existence  of  a  low-momentum  region  between  these 
vortex  paths.  Recalling  that  the  spatial  signature  of  hairpin  vortex  packets  in  streamwise- 
spanwise  measurement  planes  yields  a  qualitatively  similar  effect,  particularly  the  existence  of 
spanwise-separated  counter-rotating  wall-normal  vortices  between  which  a  region  of  streamwise 
momentum  deficit  exists,  these  patterns  in  ensemble- averaged  XciUy/\u}y\  may  be  indicative  of 
a  preferential  alignment  of  such  structures  due  to  the  roughness  below  the  measurement  plane 
or  the  generation  of  such  structures  directly  by  the  roughness.  Regardless  of  the  origin,  this 
evidence  supports  that  vortical  structures  are  likely  responsible  for  the  streamwise-elongated 
character  of  the  mean  velocity  defect  as  well  as  certain  components  of  the  Reynolds  stresses 
deep  within  the  roughness  sublayer. 
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Figure  25:  Contour  maps  of  ensemble- averaged  wall-normal  Reynolds  normal  stress,  (u/2)*h,  at 
y  =  0.0476.  (a)  Smooth;  (b)  M  =  5  model;  (c)  M  =  16  model;  (d)  Full  surface.  Green  background 
contours  highlight  local  roughness  topography  beneath  the  wall-parallel  measurement  plane. 
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Figure  26:  Contour  maps  of  ensemble-averaged  spanwise  Reynolds  normal  stress,  (u>'2)+,  at 
y  =  0.0476.  (a)  Smooth;  (b)  M  =  5  model;  (c)  M  =  16  model;  (d)  Full  surface.  Green  background 
contours  highlight  local  roughness  topography  beneath  the  wall-parallel  measurement  plane. 
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Figure  27:  Contour  maps  of  ensemble-averaged  Reynolds  shear  stress  component  ( at 
y  =  0.047<5.  (a)  Smooth;  (b)  M  =  5  model;  (c)  M  =  16  model;  (d)  Full  surface.  Green  background 
contours  highlight  local  roughness  topography  beneath  the  wall-parallel  measurement  plane. 
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Figure  28:  Contour  maps  of  ensemble-averaged  Reynolds  shear  stress  component  (tz'u/)  at 
y  =  0.047(5.  (a)  Smooth;  (b)  M  =  5  model;  (c)  M  =  16  model;  (d)  Full  surface.  Green  background 
contours  highlight  local  roughness  topography  beneath  the  wall-parallel  measurement  plane. 
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Figure  29:  Contour  maps  of  ensemble- averaged  contributions  of  ejections  and  sweeps  to  the 
mean  Reynolds  shear  stress  with  a  threshold  of  H  =  4  at  y  =  0.047J.  (a.b)  M  =  5  model; 
(c,d)  M  =  16  model;  (e,f)  Full  surface. 


5.5  Two-point  velocity  correlation  coefficients 

Finally,  two-point  correlation  coefficients  of  the  three  velocity  components  resolved  in  the 
streamwise-spanwise  measurement  plane  are  used  to  assess  the  average  spatial  structure  of 
the  flows  in  the  y  =  0.047<S  measurement  plane  within  the  roughness  sublayer.  Note  that  these 
correlations  are  computed  assuming  homogeneity  in  the  streamwise  direction  but  inhomogeneity 
in  the  spanwise  (z)  direction  to  assess  the  variations  of  these  correlations  on  spanwise  position. 
Figures  31(a)  and  31(b)  presents  puu  for  the  M  =  16  model  and  the  full  surface  flows,  respec¬ 
tively,  for  a  spanwise  position  coincident  with  the  low-momentum  region  identified  in  the  mean 
velocity  defect  for  these  two  cases.  Both  results  are  quite  elongated  in  the  streamwise  direction 
(the  puu  =  0.3  contour  extends  approximately  <5,  for  example)  and  a  region  of  weaker,  negative 
correlation  is  evident  outboard  of  the  positive  correlation  region  in  the  spanwise  direction.  The 
streamwise-elongated  nature  of  puu  in  smooth-wall  turbulence  has  been  previously  linked  to 
the  low-momentum  regions  generated  by  hairpin  vortex  packets  and  the  negative  correlation 
outboard  of  the  positive  region  in  the  spanwise  direction  represent  the  imprint  of  spanwise- 
alternating  low-  and  high-momentum  regions.  [47,  48]  In  addition,  Wu  and  Christensen  [48] 
recently  reported  similar  characteristics  of  puu  for  flowT  over  irregular  roughness  near  the  outer 
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Figure  30:  Contour  maps  of  ensemble- averaged  signed  wall- normal  swirling  strength,  (A^)  nor¬ 
malized  by  uT  and  5  at  y  =  0.0475.  (a)  M  =  5  model;  (b)  M  =  16  model;  (c)  Full  surface.  Green 
background  contours  highlight  local  roughness  topography  beneath  the  wall-parallel  measure¬ 
ment  plane. 


edge  of  the  roughness  sublayer,  supporting  the  relative  immunity  of  these  vortical  structures 
as  well  as  their  coherent  organization  into  larger-scale  structural  entities  to  roughness.  The 
present  results  lend  additional  support  to  the  contention  that  hairpin  vortex  packets  exist  in 
rough-wall  turbulence  as  well,  including  at  the  wall-normal  position  of  the  present  measure¬ 
ments  ( y  =  0.0475)  which  is  deep  within  the  roughness  sublayer  where  roughness  effects  are 
most  dramatic. 

Interestingly,  the  characteristic  streamwise  and  spanwise  spatial  scales  of  the  M  =  16  puu 
match  those  of  the  full  surface  extremely  well.  The  only  notable  difference  lies  in  the  regions 
of  negative  correlation  which  are  more  streamwise-elongated  in  the  full-surface  flow  than  in  the 
M  =  16  result.  Nevertheless,  the  consistency  between  the  M  —  16  and  full-surface  results 
highlight  the  ability  of  the  M  =  16  model  to  reproduce  the  average  structural  characteristics 
of  the  full-surface  flow  despite  not  containing  the  smaller-scale  topographical  details  of  the  full 
surface.  For  reference,  the  M  =  5  puu  computed  on  the  low-momentum  region  identified  in  its 
mean  velocity  defect  at  z  ^  0.35  shows  similar  streamwise  elongation  as  well  as  spanwise-offset 
regions  of  negative  correlation,  but  it  is  not  included  here  since  this  spanwise  position  does  not 
coincide  with  that  of  the  M  =  16  and  full  surface  low-momentum  regions. 

Figures  32(a)  and  32(b)  present  two-point  correlation  coefficients  of  wall-normal  velocity. 
pvv,  for  the  M  =  16  model  and  the  full  surface,  respectively,  at  the  same  spanwise  location 
as  puu  (i.e.,  coincident  with  the  spanwise  position  of  the  low-momentum  regions  identified  in 
their  respective  mean  velocity  defects).  Previous  studies  have  shown  that  pvv  is  much  smaller 
in  spatial  scale  than  puu  and  is  thus  more  consistent  with  the  imprint  of  the  individual  vortices 
present  in  the  flow  than  with  the  streamwise  alignment  of  such  structures  into  larger-scale 
packets.  [47,  48]  In  addition,  Wu  and  Christensen  [48]  recently  reported  the  smaller  spatial 
scales  of  the  flow  to  be  quite  insensitive  to  roughness  effects  as  reflected  in  consistency  between 
smooth-  and  rough- wall  pvv  near  the  outer  edge  of  the  roughness  sublayer.  A  similar  consistency 
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Figure  31:  Two-point  correlation  coefficients  of  streamwise  velocity,  puu ,  in  the  x  —  z  plane  at 
y  =  0.U475  at  the  spanwise  position  of  the  low-momentum  region  noted  in  the  mean  velocity 
defect  (figure  22).  (a)  M  =  16  model;  (b)  Full  surface. 


is  noted  herein  between  the  M  =  16  and  full  surface  pvv  at  y  =  0.047(5  as  both  are  quite  compact 
in  the  spanwise  direction  and  show  only  a  slight  elongation  in  the  streamwise  direction,  though 
nearly  an  order  of  magnitude  smaller  than  puu.  Thus,  the  M  =  16  model  is  able  to  reproduce 
the  average  spatial  characteristics  of  the  vortical  structures  that  populate  the  full-surface  flow 
within  the  roughness  sublayer  quite  well. 

Figures  33(a)  and  33(b)  present  two-point  correlation  coefficients  of  spanwise  velocity,  pww, 
for  the  M  =  16  model  and  the  full  surface,  respectively,  again  at  the  same  spanwise  location  as 
puu  and  pvv.  A  slight  enhancement  in  the  streamwise  coherence  of  pww  is  noted  in  the  M  =  16 
result  compared  to  that  of  the  full  surface.  Nevertheless,  the  general  spatial  characteristics  of 
the  full-surface  pww  are  clearly  reproduced  by  the  M  =  16  surface. 

Figures  34(a)  and  34(b)  present  the  negative  of  two-point  cross-correlation  coefficients  of 
streamwise  and  wall-normal  velocity,  —pUv,  for  the  M  =  16  model  and  the  full  surface,  respec¬ 
tively.  Previous  studies  have  found  puv  to  reflect  both  the  smaller-scale  nature  of  pvv,  likely 
the  imprint  of  individual  Reynolds-shear-stress-producing  events,  as  well  as  the  larger-scale 
streamwise  coherence  of  puu  indicative  of  the  collective  induction  of  such  events  by  the  multiple 
vortices  within  hairpin  vortex  packets.  [48]  The  present  results  are  consistent  with  this  behavior 
as  —  puv  is  streamwise-elon gated  and  contains  spanwise-offset  regions  of  negative  correlation. 
The  alternating  sign  of  this  correlation  in  the  spanwise  direction  is  regarded  as  the  imprint 
of  spanwisc-alternating  low-  and  high-momentum  regions,  the  former  of  which  have  been  pre¬ 
viously  found  to  embody  strong  ejection  events  while  the  latter  embody  strong  sweep  events. 
When  comparing  the  M  =  16  and  full  surface  results,  strong  agreement  is  noted  both  in  stream- 
wise  and  spanwise  extent  as  well  as  in  the  magnitude  of  the  correlation  both  in  the  positive  and 
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Figure  32:  Two-point  correlation  coefficients  of  wall- normal  velocity,  pvv ,  in  the  x  —  z  plane  at 
y  =  0.047^  at  the  spanwise  position  of  the  low-moment  uni  region  noted  in  the  mean  velocity 
defect  (figure  22).  (a)  M  =  16  model;  (b)  Full  surface. 


negative  correlation  regions.  Thus,  the  M  =  16  model  reproduces  the  average  spatial  character 
of  the  dominant  Reynolds-shear-stress-producing  events  of  the  full-surface  flow  quite  well. 

Finally,  it  should  be  noted  that  these  correlation  coefficients  computed  at  other  spanwise 
locations,  including  at  the  edges  of  the  low-momentum  regions  identified  in  the  mean  velocity 
defect  results  as  well  as  in  identified  regions  of  high  momentum,  display  similar  consistency 
between  the  M  —  16  model  and  the  full  surface.  As  such,  these  results  highlight  the  efficacy  of 
the  M  =  16  model,  containing  the  larger-  and  intermediate-scales  of  the  full- surface  topography, 
in  reproducing  the  spatial  structure  of  the  full-surface  flow  deep  within  the  roughness  sublayer 
where  the  impact  of  roughness  is  the  most  obvious.  While  not  shown,  the  M  =  5  correlations 
show  good  qualitative  consistency  with  the  M  =  16  and  full-surface  results;  however,  some, 
albeit  somewhat  weak,  quantitative  differences  are  notable,  indicating  that  the  intermediate 
topographical  scales  not  captured  by  the  M  =  5  model  have  a  measurable  impact  on  the  flow7 
within  the  roughness  sublayer. 

6  Summary 

Low-order  representations  of  highly-irregular  surface  roughness  replicated  from  a  damaged 
turbine  blade  wrere  tested  in  a  zero-pressure-gradient  turbulent  boundary  layer  under  both 
developing-  and  developed  flow  scenarios.  For  the  case  of  developing  flow,  a  model  embodying 
the  first  16  of  383  modes  (4.2%  of  the  total  modes)  which  captures  95%  of  the  full  surface  content 
accurately  reproduces  the  flow7  characteristics  of  flow  over  the  full  surface,  including  the  mean 
velocity  profile,  Reynolds  normal  and  shear  stresses,  pdfs  of  RSS-producing  events  and  quad- 
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Figure  33:  Two-point  correlation  coefficients  of  spanwise  velocity.  pWu>,  in  the  x  —  z  plane  at 
y  =  0.0475  at  the  spanwise  position  of  the  low-momentum  region  noted  in  the  mean  velocity 
defect  (figure  22).  (a)  M  =  16  model;  (b)  Full  surface. 


rant  contributions  to  the  mean  RSS.  In  contrast,  a  5-mode  surface  model  fails  to  reproduce  the 
flow  characteristics  of  the  full  surface  despite  containing  71%  of  the  full-surface  topographical 
content  and  having  an  RMS  roughness  height  that  is  84%  that  of  the  full  surface. 

For  the  case  of  developed  flow,  wherein  the  internal  roughness  layer  has  grown  to  engulf 
the  entire  bound  ary- layer  thickness  and  a  self-similar  state  is  attained  in  all  surface  cases, 
both  the  5-  and  16-mode  surface  models  faithfully  reproduce  the  character  of  the  full-surface 
flow  outside  the  roughness  sublayer,  including  both  single-  and  multi-point  statistics  as  well 
as  important  details  of  the  instantaneous  RSS-producing  events  as  discerned  from  quadrant 
analysis.  However,  this  consistency  is  a  bit  misleading  due  to  the  manner  in  which  the  flow 
statistics  are  scaled  and  compared  (by  uT  and  5).  In  fact,  since  these  flows  have  attained  a  self- 
similar  state,  this  collapse  is  simply  attributable  to  the  existence  of  outer-layer  similarity  for  the 
present  surfaces  in  accordance  with  Townsend’s  wall  similarity  hypothesis.  In  fairness,  neither 
the  5-  nor  16-mode  model  fully- reproduces  the  bulk  flow  characteristics  of  the  full-surface  flow  as 
represented  by  the  downward  shift  in  the  mean  velocity  profile  via  the  roughness  function,  A C7+, 
though  the  16-mode  model  is  quite  close  (A C/+  =  7.9  compared  to  8.2  for  the  full-surface  flow). 
Nevertheless,  this  observation  of  outer-layer  similarity  for  the  two  low-order  representations,  the 
full  surface  and  the  smooth  wall  is  actually  quite  important  because  it  highlights  the  relatively 
weak  impact  that  the  intermediate-  and  finer-scale  topographical  details  of  the  roughness  have 
on  the  outer-layer  flow.  As  such,  for  the  topography  under  study  herein,  knowledge  of  simply 
uT  and  S  would  provide  accurate  prediction  of  the  outer-layer  behavior  of  the  full-surface  and 
low-order-model  flows  from  smooth-wall  statistics. 

In  contrast,  important  differences  are  noted  between  flowr  over  the  5-  and  16- mode  models 
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Figure  34:  Two-point  cross-correlation  coefficients  of  streamwise  and  wall-normal  velocity,  — puv , 
in  the  x  —  z  plane  at  y  =  0.047(5  at  the  spanwise  position  of  the  low- momentum  region  noted  in 
the  mean  velocity  defect  (figure  22).  (a)  M  =  16  model;  (b)  Full  surface. 


and  the  full  surface  within  the  roughness  sublayer  where  the  details  of  the  roughness  are  expected 
to  play  a  pivotal  role  in  the  flow  development.  In  particular,  both  the  5-  and  16- mode  models 
yield  lower  values  of  (it'2)  in  the  roughness  sublayer  compared  to  the  full-surface  case.  Similarly, 
these  models  fail  to  reproduce  the  contributions  of  intense  ejection  and  sweep  events  noted 
in  the  roughness  sublayer  of  the  full-surface  flow.  Further,  both  low-order  surface  representa¬ 
tions  produce  an  enhanced  shortening  in  the  streamwise  extent  of  the  two-point  correlation  of 
streamwise  velocity  compared  to  the  full-surface  flow.  Taken  together,  these  observations  from 
the  streamwise  wall- normal  plane  PIV  measurements  indicate  that  finer-scale  topographical 
scales  of  the  roughness  play  subtle,  but  measurable,  roles  in  the  overall  flow  development  within 
the  roughness  sublayer  for  the  topography  considered.  In  this  regard,  it  was  noted  that  these 
finer-scale  topographical  details  captured  in  the  truncated  higher-order  SVD  modes  embody 
important  details  regarding  the  jagged  edges  of  the  large-scale  protrusions.  Thus,  the  observed 
differences  within  the  roughness  sublayer  for  the  two  surface  models  compared  to  that  of  the 
full-surface  flow  reflect  the  importance  of  these  topographical  details  in  the  near-surface  flow, 
possibly  through  modifications  of  the  separation  behavior  downstream  of  large-scale  protrusions 
and/or  vortex  shedding  from  these  features. 

The  measurements  in  a  streamwise-spanwise  plane  deep  within  the  roughness  sublayer  at 
y  =  0.047<5  highlight  how  flow  within  the  roughness  sublayer  for  all  three  rough-wall  cases  clearly 
depends  upon  the  local  roughness  topography.  In  particular,  the  full  surface  and  the  M  =  16 
model,  which  closely  resembles  the  full-surface  topography,  both  produce  similar  regions  of 
streamwise  momentum  deficit  that  are  bounded  by  intense  regions  of  (it'2),  (itV)  and  (ufwf). 
These  patterns  indicate  that  the  local  roughness  features  can  generate  significant  heterogeneities 
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in  the  turbulence  statistics  that  might  reflect  preferential  regions  of  enhanced  turbulence  pro¬ 
duction  and  possibly  dissipation.  Understanding  such  processes  is  crucial  for  the  development 
of  turbulence  models  that  can  properly  capture  the  local  impact  of  surface  roughness  on  energy 
generation,  redistribution  and  dissipation 

Finally,  the  efficacy  of  the  M  =  16  surface  model  in  reproducing  many  of  the  characteristics 
of  flow  over  the  full  surface  considered  herein  is  quite  interesting  and  provides  some  insight 
into  the  roughness  characteristics  that  have  the  greatest  impact  on  the  flow.  However,  this 
particular  level  of  modal  content  cannot  be  viewed  as  a  universal  requirement  for  constructing 
models  of  other  irregular  rough  surfaces.  Of  interest,  though,  is  the  fact  that  the  fractional 
surface  content  of  95%  for  the  M  =  16  surface  model  in  this  effort  is  consistent  with  that 
noted  in  the  recent  study  of  Johnson  and  Christensen  [33]  wherein  a  very  different  realistic 
rough  surface  was  reconstructed  with  95%  of  the  fractional  surface  content  using  the  first  20 
SVD  modes  of  the  topographical  decomposition.  Consistency  was  noted  with  flow  over  the  full 
surface  under  developing-flow  conditions  except  in  the  very  near-wall  region  Thus,  while  the 
number  of  SVD  basis  functions  included  in  a  topographical  model  will  most  certainly  vary  from 
surface  to  surface,  these  two  efforts  highlight  the  possibility  that  fractional  surface  content  might 
provide  at  least  partial  guidance  in  determining  the  appropriate  modal  content  for  low-order 
models  of  irregular  roughness.  Further  study  of  a  range  of  realistic  rough  surfaces  is  needed  to 
critically  assess  this  possibility. 
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